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A self-contained pedagogical introduction to the functional Schrodinger picture method of many- 
body theory is given at a level suitable for graduate students and also for many-body physicists 



^ , applied to various condensed matter problems including the electron gas, the Hubbard model, the 

• BCS superconductivity, and dilute bose gas within the variational approximation. It is shown that 

I the variational approximation with the Gaussian trial functional invariably leads into the Hartree- 

Fock results for both zero and finite temperature cases. 

' In order to go beyond the Gaussian results, concepts of variational and optimized perturbation 

, theories are introduced using simple quantum mechanical languages. Then, the variational pertur- 

__ ■ bation theory in the functional Schrodinger picture is applied to the A(^'* model to yield the effective 

'"7^ ' potential to the third order. The optimized perturbation theory is also applied to the A(/!>* model to 

I , yield a general expression up to the second order. 
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FOREWORD 



These notes represent the lectures given at Johannes Kepler Universitat-Linz, Austria in October and November, 
2002 as " Special Topics in Many-Body Theory" . The intent of these lectures was to provide a self-contained tutorial 
introduction to the functional Schrodinger picture approach to many-particle systems. 

Condensed matter systems are basically many-particle systems. Whether the system in consideration is a crys- 
talline solid, an amorphous material, or a liquid, one has to deal with substantially large number of particle which 
often amounts to 10^'^/cm'^. This large number still eludes exact numerical calculation even with ever increasing 
computational capability of modern day computers. 

Many-body theories based on non-relativistic quantum field theories are conceived to deal with almost infinite degrees 
of freedoms of many-particles in condensed matter systems. Many-body theory closely followed its relativistic brother, 
quantum field theory in its development. In quantum field theory, there exist basically three alternative approac;hcs, 
namely operator (Heisenberg) , functional integral (path integral), and Schodinger approaches. Comparing to the 
former two methods, the last one Schodinger picture approach was early regarded less convenient and clumsy for 
development of field theories. However, in 80's and 90's, theorists have found that it offers some distinct advantages 
over others when suitable schemes are employed. Since then the Schodinger picture has been drawing wide interests 
and field theories based on the functional Schodinger picture have been developing rapidly. 

In condensed matter physics, so far the Schodinger picture scheme has been much neglected. Only very recently, 
efforts to develop the Schodinger picture method into a useful many-body theoretical tool has started. In this lecture, 
we present such an effort and try to give a pedagogical overview on this method. 

Since these notes are designed for beginning graduate students and newcomers in this field, we tried to make the 
notes selfcontained. In this respect, we tried to limit the number of references to those essential to read these notes. 
Other relevant references can be found in the references of the papers and books cited in these notes. Although most 
of essential steps are described in detail to help the readers, some of the derivations are left as problems, so that 
readers can gain more insights and expertise through exercise. 

I. INTRODUCTION 

Condensed matter systems are basically many-particle systems with many particles (~ 10^^ per cm'^). Now, the 
basic question is whether it is possible to solve the many-particle Schrodinger equation directly? 

JV 1 ''^ 

H = j2T{xk) + ^ E nxk,xi), (1) 

k=i fe/;=i 

d 

■■■XN,t)= H'^ixi ■■■xn, t). (2) 

If we have 4'(.ei • • ■ XN,t), then we can calculate all the physical quantities we need. 

Expansion of ^'(a;i • • ■ XN,t) can be written in terms of the single-particle wave function {xk) 

^{x^---XN,t)= C{E[---E'^,t)^E[{xi)---^E'JxN), (3) 

where E[ ■ ■ ■ E'^^ represc^nt single particle quantum numbers. C{E[ ■ ■ ■ E'j^,,t) is the coefficient for a particular c;onfig- 
uration. In general, it is not practical to continue to use ■ ■ ■ XN,t) with N ~ 10^'^. Therefore, a suitable scheme 
to handle the large degree of freedom has to be adopted. Here, we resort to the technique of second quantization. 

A. Second Quantization 

There are several ways for introducing second quantization to the many-body Schrodinger equation. An explicit 
and detailed treatise is given in Chapter 1 of Fetter & Walecka (Gen. Ref. 3). However, here, we follow a rather 
concise but less detailed path to the second quantization commonly adopted by field theory textbooks (Chapt. 2, 
Gen. Ref. 1 & §12 & 13, Gen. Ref. 4). 

The Schrodinger Wave Field of Bose Statistics 
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(1st quantized) Schrodinger equation is given by 



The conjugate complex equation is 



V^* + = i/i*. (4) 

2m 



V^** + V(x)^* = -im*. (5) 

2m 



Find a Lagrangian such that the corresponding Euler-Lagrange equation leads back to Schrodinger equation. 

d fdL\ dL 



, 1 . , i = 1 • • • / (6) 



Such a Lagrangian is given by 

l;2 



which leads directly to 

d 5L SL r_. 



- V{x)^ + i^V^^ } = 0. (8) 



dtS-^* S^* { ^ ' 2m 

(Note: The above Lagrangian is not Hermitian. We can work with a Hermitian Lagrangian (See Chap. 2, Gen. Ref. 

1.) at the expense of slight complications.) 

Define the canonically conjugate momentum in the usual way, 

7r= -J- 9 

The Hamiltonian is given by 

H = j{TT^- C)d^x 

= y"**(f)|-^V2 + F(x)|*(f)d3x. (10) 

In the above Hamiltonian, only wave nature appears. In order to restore the particle nature, we expand ^ in terms 
of the 1st quantized eigenfunctions, 

2m ) 



The expansions are 



i)^{x,t) = e-^^-'xl^^ix). (11) 
*(f) = 5^6^(t)V^(f), 

^*{x) = J2bUt)r^{^), (12) 



where b^{t) = 6/,(0)e"S"^''*. 
With these expressions, we obtain 
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H = Y,E^b%. (13) 

In order to quantize the wave field \E', we require 

[7r(f),*(f )] = ). (14) 

(Remember the 1st quantization postulate \p,x] = j-) 
Thus, we obtain 

\^ix'),¥{x)]=5{x-x'). (15) 
Also we have additional commutation relations 

[*1'(x),*^(:r')] =0. (16) 

We find that and fej^ satisfy 

[^,&t]=V. K,b.]=0, [blbl] = 0. (17) 
With the new Hamiltonian operator, the 'Schrodinger equation' is given by 

= E<i>, 

^Y,E^b%, (18) 



t,2 

-—V^ + V{x)\'^{x)cfix 



where 

^=n^(^D"''*0' (19) 

E = ^E^n^ = 0,1,2, •• •. 

A series of energy levels E^ are occupied by a definite number of quanta. This shows clearly that field quantization 
guarantees the corpuscular nature of the Schrodinger wave field. 

Expectation Values in the New Formalism 

(How do we relate the new theory to reality?) 

► The particle density operator 

**(a;,f)*(a;,t) — > *^(a;, t)«'(a;, i) (20) 
(1st quantized) (2nd quantized) 

Pop{x) = ^\x)'^{x) : particle density operator 

p{x) = m^\x)^{x)\^) (21) 

<Problem> Show that the above p{x) gives an identical result to the 1st quantized particle density p{x) = ip*{x)tp{x), 
when only a single particle is present in state k. 

t The position operator 

In the 1st quantized form, we have 
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a; = y il)*{x,t)x'tjj{x,t)d^x. 

The second quantized form is given by 

Xop = j ip^x)xil){x)d^x, (22) 

X= {^\Xop\^) 

= ($1 j ■tlj^{x)xi){x)d^x\^). (23) 

• The potential energy operator 

Vop = J ij\x)V{x)i;{x)d^x. (24) 

• The interaction energy operator 

V{x,x') = V{x-x'). 

In the first quantized form, it is given by 

^ J p{x, t)V{x, x')p{x', t)d^xd^x' 
= ^*{x,t)^{x,t)V{x - x')^*{x' ,t)^{x' ,t)d^xd^x' . (25) 

In the second quantized form, the order of the operators is important. 

When there is one or zero particle present, the interaction energy is zero. 
Thus, we should have \E' on the right and on the left, 

^ j 'i/^x)'^'^ {x')V{x - x')'^{x')'^{x)d^xd^x' . (26) 
The expectation value is given by 

j ^^/\x)'^^{x')V{x-x')'^{x')'^{x)d^xd^x'\<^>). (27) 

<Problem> Calculate the expectation value of the Coulomb interaction energy V{x — x') = e'^/\x — x'\ for $ = bj^^^o 
and$ = 6j,^6t^$o • 

Quantization of the Schrodinger Wave Field of Fermi-Dirac Statistics; Fermions 

In reality, electrons which are main constituents of the many-body problems in condensed matter are fermions, 
which obey F-D statistics. 



We expand the wave field into 



vl/t(a;)=^at^;(x), (28) 



where i)fj,{x) satisfies 



-^V^+V{x)]^l;^{x)=EM^). 
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In order to satisfy the F-D statistics (Pauli exclusion principle), they should obey the following commutation relations. 

{ai,at} = alal + alal = 

{at,a.} = V (29) 
{an,a„} = 

Using these relations, we obtain 

{^^x),^ix')} = S{x-x'), 

{*t(^),vi,t(^')} = 0, (30) 
{*(a;),*(x')} = 0. 

<Problem> Derive the above commutation relations using the commutation relations of aj^{x) and ai,{x). 
The Hamiltonian operator is given by 

H = J ^^x)l^-^V^ + V{x)^^{x)d^x (31) 
The corresponding Schrodinger equation is 

= E^. (32) 



We obtain as eigenfunctions 



with N = J:^ N,, and E = J2^ E^n^. 
First consider a one-particle state 



= Wialr-^o, (33) 



$ = ^c^aj,$o (34) 

which is a superposition of one-particle states ajj^'o- 
We replace aj, by \E''^ using 

j V^(a;)*t(a;)d3 X = a' . (35) 

/ ^c^^/;^(x)*^(a;)d3a;$o 

= j f{x)¥{x)d^x^o (36) 

$ represents a superposition of the single particle, which maybe thought as having been created at different places 
with a probability amplitude f{x). 

<Problem> Show that '^^{x') creates a particle at the place x'. 
<Sol.> ^^x)^{x)^^x')^o = ^^x)6{x - x')^o = S{x - a;')*^(a;')*o- 

We now have 
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= J J (-^V^ + V{x)^ f{x')S{x - x')d^xd^x'^o 



) ) fix)d'x. (37) 
Here, we use the relation = = E J f{x)^^x)^Qd^x. This is valid only when 

.^V'+Vix)^f{x) = Ef{x). (38) 

This one-particle example shows that the second quantization formalism is identical to the usual Schrodinger theory. 
• What about a 2-particle state? 
Consider the most general two-particle state 

^=Y1 c,*i,*2 41^2^0, (39) 

which can be rearranged 

$= / ^{J2i^^^Ax)Mx')c^.^.2)^Hx)'^Hx')^od^xd'x' 

f{x,x')^\x)^'^{x')^od^xd^x' . (40) 
Now we show that f{x,x') is antisymmetric. Interchanging x with x', we have 

f{x',x)^^{x')<i/^x)^od^x'd^x 
= -j j f{x\x)-^\x)-^\x')^Qd'^x'd^x. (41) 
Thus, we note that f{x,x') = —f{x',x) by comparing to Eq.(40). 

Now, we assume the Hamiltonian operator of interacting electrons in second quantization is given by 

72 I 1 ,Tf/^\j3„ 



H = J *t(x) (^-^V^ + V{x)j ^{x)d-'x 

+ - [ [ ^Ux)^Ux') , „ ^(a:0^(x)d3xdV. (42) 
2 J J \x-x'\ 

Combining Eq.(40) with (42), we obtain 

<Problem> Prove Eq.(43). 
► n-particle state 

If we apply the field Hamiltonian operator to an n-particle state, 

^ = j d^xi ■ ■ ■ d^Xnfn{xi,-- ■,Xn)^\xi) ■ ■ ■ *1'(a;„)$o, (44) 



we will find the n-body wavefunction must satisfy 



J I J 



In this way, aU of the n-body first quantized systems can be contained in the corresponding second quantized theory 
(quantum field tlieory). In this sense, 'many-body theory' is basically a 'nonrelativistic quantum field theory'. 
For interacting electron systems, the Hamiltonian is given by 

+ j d^xd^x''i>\x)^\x')V{x,x')'^{x')'^{x) (46) 
= ^al{r\T\s)as + ^ J2 alal{rs\V\tu)auat. (47) 

r,s r.s,t.u 

(Note the order of operators.) 



B. Pictures in Quantum Mechanics 



Schrodinger Picture 

The usual elementary description of quantum mechanics assumes that the state vectors are time dependent, whereas 
the operators are time independent. The Schrodinger equation in the Schrodinger picture takes the form 

in^^\^s{t))=H\^s{t)). (48) 

A formal solution is given by 

|*5(<)) = e-'«(*-*«)/''|M'5(to)). (49) 

Generally, undergraduate quantum mechanics dwells on the Schrodinger picture formalism, although the Heisenberg 
picture, which we will discuss shortly, is regarded as a convenient tool for a formal development of quantum mechanics. 
However, in many-body physics, the situation is completely reversed. It will be shown that the Heisenberg and the 
interaction pictures naturally lead into perturbative expansions of the physical quantities, which are the basis of the 
Green's function formalism. 



Therefore, almost all the many-body theoretical schemes are based on the Heisenberg (interaction) picture. Only 
recently, several quantum field theorists found that the Schrodinger picture provides a convenient scheme to calculate 
the physical quantities in selected quantum field theoretical problems. Last few years, the authors and their collab- 
orators started to develop a many-body theory based on the Schrodinger picture and have succeeded to provide a 
general theoretical framework. In this lecture, the many-body theoretical scheme based on the functional Schrodinger 
picture will be introduced. 



Heisenberg Picture 



In the Heisenberg picture, we pass the burden of the time-dependency from the wavefunction to the operator. The 
state vector (wavefunction) in the Heisenberg picture is defined as 

\^H{t)) = e'^*/'^|*s(t)). (50) 

Gombining Eq.(49) with to = 0, we obtain 

|*ff(i)> = l*5(t = 0)) (51) 

which is independent of t. 

Gonsidering the fact that the expectation value should remain same in both pictures. 
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{^s{t)\Os\^s(t)) 

= (*h|Oh|*h), (52) 

we obtain 

Oh (i) = e'"*'^Ose-'"^'^ . (53) 

We observe that the state vector in the Heisenberg picture is now independent on t, whereas the time-dependency is 
now carried by the operator. This fact is summarized as follows 



Schrodinger picture Heisenberg picture 
SteteKrt T-dep. T-ind. 

Operator T-iud. T-(k>l). 

The time derivative of Eq. (53) yields the Heisenberg equation of motion, 



in^OH{t) = [OH{t),H\. (54) 



<Problem> Prove Eq.(54). 



This result determines the equation of motion for any operators in the Heisenberg picture. In particular, if Og 
commutes with H , then Oh is a constant of the motion. 

Interaction Picture 

When the Hamiltonian can be expressed as 

H = Ho + Hi, (55) 

where Hq yields a soluble problem, we can introduce the interaction picture, which is convenient for perturbative 

treatment of many-particle systems. 

Define the interaction state vector in the following way 

\^i{t)) = e'"°'/''\^s{t)) (56) 

(* Note that in the exponent, it is not H but .ffo-)- Then 

We, therefore, obtain the following equations. 

zn^^\^j{t))=Hj{t)\^j{t)) (57) 
Hi{t) = e'"ot/njj^^-iHot/h ^58^ 



For any operators. 



Thus, we have 



{^s{t)\Os\'fs{t)) 

= (*/(t)|e^^«*/'*Ose-'^°*/'^|*/(i)). (59) 



Oi{t) = e»^^ot/R(9^e-ii?ot/ft_ (60) 
Differentiating the above equation with respect to time, one obtains 
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= [Oi{t),Ho]. (61) 

We now solve the equation of motion in the interaction picture. Define a unitary operator that determines the state 
vector at time t in terms of the state vector at the time 

\^i{t)) = U{t,to)\^j{to)). (62) 

U satisfies 

U{to,to) = l. (63) 

Using the Schrodinger picture Eq.(49), 

\^/j{t)) = e*^»*/'*|4's(f)) = e'^°'/'^'e-*-f^(*-*°)/'*|*s(to)) 

= e'^''*/Re-'^(*-*'')/'^e-'^°*"/'*|*/(io)), (64) 

which identifies 

U{t, to) = e*-f^o*/'ie-*^(*-*«)/''e-*-f^"*»/''. (65) 

We can immediately show that 

1. W{t, to)U{t, to) = U{t, to)W{t, to) = 1 
which implies U is unitary 

UHt,to) = U-\t,to). (66) 

2. U{tut2)U{t2M) = U{tutz). (67) 

3. U{t,to)U{to,t) = l. (68) 

Although Eq.(65) is the formal solution of Eq.(62), it is not very useful. 
Prom Eq.(57), one can write 

inyj{tM)\'^i{h)) = Hi{t)uit,to)\'^i{to)), 

which is true for arbitrary |^'j(fo)). 
Thus, we obtain 

ih^^U{t, to) = Hi{t)U{t, to). (69) 
Integrating this equation from to to t, it is shown 

U{t,to) - U{to,to) = f dt'Hi{t')U{t',to), 

" J to 

which gives the relation, 

U{t, to) = l-- dt'Hi{t')U{t', to). (70) 

" -Jto 

Thus, we arrived at an integral equation for the unitary operator U(t,to)- 
The above equation is solved by iteration. 

U{t,to) = 1 + (^^) j^* dt'Hiit') + (^^^ ^* dt' ^* dt"Hj{t')Hi{t") + ■■■ (71) 
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Consider the third term in this expansion. It may be written as 

ft nt' 



I dt' I dt"Hi{t')Hi{t") 

J to J to 



\ f dt' f dt"Hi{t')Hi{t") + \ f dt" [' dt'Hj{t')Hi{t") 

^ J to Jtn ^ Jto Jt" 



(72) 




Now the second term becomes (with t' ^ t' 

rt rt 

2 

Thus, we have altogether 

rt rt' 



I f dt" f dt'Hi{t')Hi{t") = ]- f dt' f dt"Hi{t")Hi(t'). 

^ Jto Jt" ^ Jto Jt' 



f dt' f dt"Hi{t')Hi{t") 

Jto Jto 



to J to 



dt' / dt"[Hi{t')Hi{t")9{t' -t")+Hi{t")Hi{t')9{t" -t')], 



(73) 



where 0{t) is the step function. We note that the operator containing the latest time stands farthest to the left. We 
call this a time-ordered product of operators, denoted by the symbol T. Thus the above expression can be rewritten 



f dt' [ dt"Hi{t')Hi{t") = ]-[ dt' j dt"T[Hi{t')Hi{t")]. 

to to J to J to 



Finally, a general expression for U is obtained 

oo / ■ \ " 1 rt 



/ _ ■ \ " 1 ft ft 



(74) 



(75) 



Now, a Heisenberg operator can be related to an interaction operator using U. 

_ ^iHL/h^-iHot/hQ^^iHot/h^-iHt/h 



Eq.(65) yields 



OH{t) = U{0,t)Oi{t)U{t,0). 



(76) 
(77) 



In addition 
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= |*s(0)) = |f /(O)) , Os = Oh{0) = Oi{0), (78) 

so that all three pictures coincide at time t = Q. 

These relations provide a formal way to obtain the exact eigenstates of the interacting system from the noninteracting 
ones. 



Since l^'j?) is independent of time, it satisfies the time- independent form of the Schrodinger equation 

H\^H)=E\<bH). (79) 

These state vectors arc therefore the exact eigenstates of the system and are naturally very complicated for an 
interacting system. This can be expressed 

\^h) = |*/(0)) = f/(0,io)|*/(io)). (80) 

If we choose [^^/(to)) = oo)) and H = Hq + e^^l*l_ff/ and use the concept of "adiabatic switching on", we can 

express an exact eigenstate of the interacting system in terms of an eigenstate of Hq. 

= |*/(0)) = U,{0, -oo)|*o), (81) 
Ho\^o) =Eo\^o)- (82) 

(For details, see Chapt. 3 of Gen. Ref. 3.) 

The above relation allows a convenient and systematic calculation of physical quantities using Green's function method 
in the Heisenberg picture. 



C. Green's Function Method 



Zero- Temperature Formalism 

The concept of a Green's function (or propagator) plays a fundamental role in the many-particle theory. 
The single particle Green's function is defined by 

zG^^ixt, xt) = j^-^ , (83) 

where |\l/o) is the Heisenberg ground state of the interacting system satisfying 

iJ|*o) = E\^q) (84) 

and tpHa{x,t) is a Heisenberg operator, 

ijHa ix,t) = e'"'/"^^ (a;)e-^^*/'^ . (85) 

Here, the indices a and /? label the components of the field operators; a and /3 may represent two values for spin ^ 
fermions, whereas there are no indices for spin zero bosons. 
The T product represents the time ordering, 

T[lPHa{xt)^%fs{x't')]= iJHc{xt)^%fs{x'if) t > t' , 

±i}^jjp{x't')i)Ha{xt) t'>t, (86) 

where the upper (lower) sign refers to bosons (fermions). 

Here, we ask why the Green's function is so useful. The answers are as follows: 



1. The Green's function offers a convenient basis for the perturbation expansion in combination with the Feynman 
rules. 
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2. Although the ground-state expectation value implies the loss of much detailed information about the ground 
state, still the single-particle Green's function contains the observable properties of great interest: 

(a) The expectation value of any single-particle operator in the ground state of the system. 

(b) The ground-state energy of the system. 

(c) The excitation spectrum of the system. 

The single-particle operator is expressed as 



where 



Jl3a{x) is the first-quantized operator. 

The ground-state expectation value of the operator density is given by 



(*o|*o) 
lim J0a{x) 

X — >X ' ^ 



= ±i lim \imy^Jaa{x)Ga3{xt,x't') 

a/3 

= ±zlim lim tr[J(a;)G(a;t,a;Y)]. 

t'^t+ x'—^x 

For example, the number density operator, n{x) = 'il)\x)'ij){x), becomes 

{n{x)) = ±itvG{xt,xt+). 

The kinetic energy (T) is 



(T) =±i [ d^x lim 

J x'^a 

and the potential energy is 

{V) =±\i I d^x lim lim V 



VhT:G(xt,x't+) 

2m V > y 



(87) 



(88) 



(89) 



(90) 



(91) 



(92) 



<Problem> Prove Eq.(92). 

The total ground state energy is solely given in terms of the single-particle Green's function. 



E={T + V) = (H) 

1 f ^ . 
i—i I d X lim lim 



ih 



dt 



2m 



tTG{xt, x't') 



(93) 



In addition to the above useful quantities, we can also obtain information on the excitation spectrum of the system 
[Section 7, Gen. Ref. 3. ]. 

Now the usefulness of the above relations boils down to the point that how accurately we can calculate the Green's 
function. In this lecture, it is not intended to explain the Green's function method. Here, we only mention that the 
Green's function method leads into a perturbative expansion without proof. 
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The exact Green's function may be written as 



/oo 
-OO 



(*o|^|*o) 



(94) 



Here, S = U^{oo, — oo). The above expression clearly shows that we can calculate the Green's function order by order 
of the perturbation, Hi , up to the desired order. This is greatly facilitated by the Feynman diagram method. 

Dyson's equation 

The concept of the Green's function method can be beautifully represented by the Dyson's equation pictorially. 

G = Go + GoSGo, (95) 

where S represents the self-energy due to the interaction and Go the free particle Green's function. This equation is 
graphically represented. 




Using the concept of the proper self-energy, this expression can be modified as 

G = Go + GoT,*G, 

which is graphically given 



+ ^ 



(96) 




A proper self-energy is the self-energy insertion that can not be separated into two pieces by cutting a single particle 
line. 

Eq.(96) and its graphic representation clearly shows the perturbative nature of the formal solution as follows 



G = Go + Go^*Go + Go^* GqYI* Go + * 



(97) 



This perturbation expansion allows one to calculate the Green's function and consequently associated physical ob- 
servables as accurately as possible when the perturbation (or the mutual interaction) is weak. 

Finite- Temperature Formalism 

At finite temperatures, it is necessary to separate the calculation into two parts. The first step is the introduc- 
tion of a temperature Green's function Q. This function has a simple perturbation expansion similar to that for G 
at T = and also enable us to evaluate the equilibrium thermodynamic properties of the system. The second step 
then relates 5 to a time-dependent Green's function that describes the linear response of the system to an external 
perturbation and provides the excitation energies of the system containing one more or one less particle. 

At finite temperature, it is most convenient to use the grand canonical ensemble, whose grand Hamiltonian is given 

by 



K = H-iiN. 



(98) 



The grand partition function is written as 
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-p(E,-tiN) 



N j 

= EEw>"'^''"'^''^i^.) 

N j 

= TV (e-^(^-''^)) 

= TV (e-^^) . (99) 

The thermodynamic potential is given by 

Q(T, V,/i) = -fcsTlnZc. (100) 

The statistical operator po is given by 

= e''(^-^). (101) 
For any operator, the ensemble average {O) is obtained with the prescription 

(O) = Tr (pgO) = Tr (e^^^^^-^)©) 

Tr (e-^^O) , , 

= (102 

We introduce the modified Heisenberg picture operator 

Ok{x, t) = e^^'^Os{x)e-^^/^. (103) 
In particular, the field operators assume the form 

V'ifa(a;,T)=e^-/Va(a;)e-^^/^ 

V-kl^^ r) = e^^/Vi(a;)e-^-/^ (104) 

Note that V'lfaC^' ''") ^'^^ adjoint of il>Ka{x,T) as long as r is real. 
The single-particle temperature Green's function is defined as 

ga0{xT,x'T') = Tr{pG7V[V'ifa(a;,T)Vlc^(a;',T')]} • (105) 
Here the symbol orders the operators with smaller r at the right. 

The temperature Green's function is useful because it enables us to calculate the thermodynamic behavior of the 
system. If H is time independent, then Q depends only on the combination t — t' and not on r and r' separately. 

By definition, 

tig{xT,XT+) = T^TT:[PG'4'Ka{x,T)lljKa{x,T)] 
a 

= Te^^^^^TV [e-''^e^-/Vi(a;)Va(a;)e-^^/'^" 

a 

= T{n{x)). (106) 

A single-particle operator can be written 

(J) = Tr(pG J) 

= / d^x lim lim Jpa{x)Ga0{xT,x' t') 

a/3 •' 

= Tld^xlim lim tr[J(a;)a(a;T,a;'r')]. (107) 
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Particular examples of interest are 



(a) =T J <fxtT[ag{xT,XT+)], 

(T) = T / d^x lim „ tTg(xT,x'T+). 
J x'^x 2m 



(108) 
(109) 



The ensemble average of two body operators usually requires a two-particle temperature Green's function. 
However, the mean potential energy of two-body interactions can be solely expressed in terms of Q as in the zero- 
temperature case. The expression is given by 

<y) = \j J rf'a;"y(a;-a;")Tr[pG^^(a;)V:l(a;")V'7(2^")V'a(a;)] 



1 f ^ 

=F- / d X lim lim 

2 J x'—^x t'^T"^ 



dr 2m 



tT:g{xT, x't'). 



(110) 



<Problem> Starting from the Heisenberg equation of motion on ipKaixr), prove Eq.(llO). 

The internal energy is given by 

E = {H) = {T + V) 



=p— / (fix lim lim 

x' — >x r' — 



dr 2m 



trg{xT, x't'). 



(Ill) 



Here, we give the expression for the thermodynamic potential in terms of the single-particle Green's function without 
proof (For details see Section 23, Gen. Ref. 3 ). 



fl{T, V, jj) = Qo{T, V,ii)t A-^dA J 



lim lim - 

ic'— r'— 2 



tvG^ixT, x't'), 



where the Hamiltonian is written as 



H{\) = Ho + XHi. 



(112) 
(113) 

(114) 



Again, here we observe that the problem of obtaining useful physical quantities boils down to the point how one can 

obtain exact Green's function. 

We end this section without going into the Green's function formalism. Instead, we just give the expression for the 
Green's function which shows the perturbative nature of the formalism. 



trg,^(xT,xV) = -Tr £ 



i-n)-^ 



n=0 



/ dn-- - / dTnTr[K^{Ti) ■ ■ ■ K^{Tn)i^io,{xT)lPip{x'T')\ 

Jo Jo 



-PKo 



E 

ra=0 



{-ny 



/ dn-- - / dTnTr[K^{Ti) ■ ■ ■ K^{Tn)\ 

Jo Jo 



(115) 



The denominator is just the perturbation expansion of e~^^. It serves to eliminate all disconnected diagrams in the 
numerator. It should be reminded that the Dyson's equations in the zero-temperature are equally valid in the finite 
temperature formalism. 

In this section, wc showed that the Green's functions in the zero and finite temperature formalism naturally lead 
into perturbative expansions. The perturbative calculations are greatly facilitated by the Feynman diagram rules. 
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As will be shown later in this lecture, in general, perturbation theory yields an asymtotic rather than convergent 
series. Under appropriate circumstances, it may yield a useful physical approximation but it cannot be applied to 
arbitrary precision. Therefore, when the perturbation, which often represents interactions between particles is large, 
then the outcome of the Green's function method is in doubt. It should be noted that many of challenging problems 
in condensed matter physics originate from strongly correlated systems. In order to handle such strongly correlated 
systems, various non-perturbative schemes have been developed. However, in many cases, the approaches lack system- 
atic means of improvement beyond the first obtained result. In later chapters, we address this problem and present 
a method, which provides a systematic improvement scheme beyond the mean field results for strongly correlated 
systems. 



II. SCHRODINGER PICTURE IN MANY-BODY THEORY 



The Schrodinger picture for many-body (field) theory is a natural extension of one- or few-body quantum me- 
chanics. 

For ordinary quantum mechanics, we start with a Hamiltonian operator and canonically quantize by postulating 
commutation relations between coordinate of position operators and their conjugate momenta. We reach the coordi- 
nate Schrodinger picture representation by representing the diagonal position operator with its eigenvalues, and use 
a differential representation of the commutators by replacing the conjugate momenta with derivatives. Coordinate 
representation of state vectors are called wavefunctions. The Schrodinger equation becomes a differential equation. 

For many-body (field) theory in the Schrodinger picture, wavefunction should carry infinite number of degrees 
of freedom. This can be accomplished by substituting 'function' by 'functional'. Differential representations of the 
canonical commutators are obtained by replacing conjugate momenta with functional derivatives. Coordinate repre- 
sentations of state vectors or elements of Fock space are wavefunctionals. The Schrodinger equation is a functional 
differential equation, whose solutions, the eigenfunctionals of the Hamiltonian functional differential operator, repre- 
sent possible states of the system. 

In addition to the functional Schrodinger picture method, the path integral formulation of many-body theory also 
requires functional calculus. In this connection, this approach is called functional integral method of many-particle 
(field) theory. In this lecture, we will not treat functional integral method. Instead, we refer to Gen. Ref. 2. 

First, we learn about functional and its calculus. 



A. Functional Calculus 



Functional 



A functional space is an infinite-dimensional space where each point in the space is a function on space-time. That 
is, each point in the function space is a mapping of space-time into the real or complex numbers. A particular point 
in a function space will map points in space-time to scalars, spinors or vectors. 

A mapping of the points in the function space to numbers is called a functional. Thus, a functional associates a 
number with each function on space-time. 

For example, let a be a point in a scalar function space A. a is a scalar function on space-time, a = a{x) G C ov TZ. 
Let F be a functional on A. F maps points in A into numbers. F = F[a\ &C ovTZ. A simple example of a functional 
is F[a] = J a{x)dx. 

Functional Derivatives 



The functional derivative is formally defined by 

Sa e^o e ^ ' 

where 6 is the Dirac 5-function distribution. Similarly, we can define a functional directional derivative in the direction 
of the function A (a;) as 

S^^^^F[a + eX]-F[a]_ 
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(ii) 




Fig. II. 1 (i) a is a function on space-time. It maps points in space-time to real or complex numbers, (ii) a is a point in 
function space A. The functional F[a] maps points in the function space A to real or complex numbers. The functional 
F[a] turns functions into numbers. 



Also by SF[a]/6a{x), we mean the change in F with respect to a change in the function a at the point x only. 
Formally 



SF[a] _ ^.^^ F[a] - F[a] 
da{x) E^o e 



(118) 



Here, a is equal to a except at the point x where d = a + e. For example, if F[a] = J a{z)S{z — y)dz = a{y), then 
F[a] = J {a{z) + e6{x - z))6{z - y)dz = a{y) + e5{x - y). 
From above, we observe that 



^ajy) ^ a{y) + e5{x - y) - a{y) 
5a{x) c^o e 
= 5{x-y). 



Similarly, 



5xa{y) 
S\a{x) 



X{x-y). 



(119) 



(120) 



Important examples include the action fimctional S = J L{qi,qi,t)dt, the Lagrangian functional L = 
J £{4>,\/(j),(j),t)d'^x and the Hamiltonian functional H = J{n{x,t) ■ dt4>{x,t) — C)d^x. In fact, in deriving Eq.(8) 
from Eq.(7) we have implicitly used the functional calculus. 

As an exercise, wc consider the action for the Ac/)"* model in the scalar field theory (We use the scalar field theory 
as a prototype to develop the functional Schrodinger theory.). The action is given by 



= j dt j d^yC{x) 



(121) 



18 



Here the notation 

ds^ = dx^dXf, = dt^ - (rfa;^ + dy'^ + dz^) 

is used with c = 1. We used 

= {x'^ , x^ , x'^ , x^) = {ct,x,y,z), 
Xfj, = {xo,xi,X2,X3) = {ct,-x,-y,-z). 

Also, the differential operators are defined 



For future reference, we also list 



D = d''d =-—-(— — — 

^ (? d1? \dx'^ dy"^ dz^ 



C2 df^ 



called the d'Alembertian operator. 

The energy-momentum 4-vector of a particle is 



given the invariant 

or, when c = 1 

Often p • a; is used for p^a;'* 

From Eq.(121), we have 



'E \ (E 
= I —,P I , Pfi= ( —,-p 



E^ 



2 i7i2 -a 2 
= E — p = m . 

p - X = p^x^ = Et — p ■ r. 



5(t>{x) 



and 



5<t){x) 



d^y\{d'^cl>d,cl>) = j dS\{ct>d^d,cl>) 

d^y6\x ~ y)d^dM = -d''d^<j){x), 



d^y^iy) =nj d^yr-Hy)6\x - y) 
= n(/)"-^(a;). 



Thus 



'^'^t'^i = - (d,d^<j^{x) + + ^ct^Hx)] = 0, 



5(t){x) yi^ J ' ^ ' 3! 
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which is the Klein- Gordon equation with the <f)^ interaction. 

Functional Differential Equations 

In the Schrodinger picture of many-body theory, the Schrodinger equation is a functional differential equation. For 
an action S[^], the field equations are 

^ - (133) 

A functional differential equation can be thought of as an infinite set of coupled partial differential equations or a 
partial differential equation on an infinite-dimensional vector space. 
For example, consider 



/*(- 



This equation is satisfied by 

F[a\ = r?exp(± j dxs/b(x)a{x)). (135) 

<Problem> Show that the above F[a\ satisfies Eq.(134). 

A more realistic form which may appear in the real problem is given by 

/ '^'^J^J^^^"^ ^ j dxdyf{x,y)a{x)a{y)F[a] = fiF[a]. (136) 

Assume f{x, y) is symmetric in x and y. 
First we substitute 

F[a] = r;exp(G[a]), (137) 

and solve the resulting equation for G[a]. 



/*(< 



5a{y)5a{y) \Sa{y) 
= 1^- j dxdyf{x,y)a{x)a{y). (138) 
To proceed, we must match the power of a{y) for both sides. Thus, we write 

G[a] = ^ j dxdyg{x,y)a{x)a{y). (139) 
G[a\ will be a solution provided the symmetric g{x,y) satisfies 

11 = j dyg{y,y) and f{x,y) = Jdzg{x,z)g{y,z). (140) 

This technique of comparing the powers of both sides is called the 'power counting' method. 

Functional Integrals - Bosonic Variables 

Functional integrals are integrals of functionals over function spaces and, in general, are not well defined. There 
are only two types that can be done exactly. The first type is the (5-functional, S[a — 



j Da5[a - ^] = 



(141) 
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D in the measure indicates that the integral is over the function space containing a{x), which can be expressed 

£)a = J|rfa(f). (142) 

Correspondingly, 

5[a-^]=^5{a{x)-m)- (143) 
Thus the integral in Eq.(141) is an infinite product of independent integrals, one for each x. Thus, we have 

j DaF[a]5[a - ^] « j J\ da{x)5{a{x) - ^(f ))-F[a] 

X 

= n / da{x)5{a{x) - (,{x))F[a] 

X 

= m- (144) 

Often it helps to think of the functions, such as a(x), as infinite dimensional (column) vectors where x plays the 
role of an index. With this view, a functional integral looks like the natural infinite-dimensional limit of ordinary 
finite-dimensional integrals. The second type of functional that we can integrate exactly is the Gaussian functional. 

j £)ae-^["l, where F[a] = j dxa'^{x) (145) 



We have 



J Dae- 1 '^'"''^''^ = J ]Jda(a;)e-Sx"'(^) 

= 11 f da{x)e-'^'^^^ =l[V^, (146) 



where we used dxe = \/^. 

J —oo V a 

Following the same steps, we have 



J Dae- 1 = H yig^. (147) 



We regard f{x) as an infinite-dimensional diagonal matrix 

f{x,y) = f{x)S{x-y). 

f{x) also represents the diagonal elements of the diagonal matrix. They are also the eigenvalues of the matrix. 
The product of the eigenvalues is the determinant of the matrix, hence 



= det-^/2 {^], (148) 



so 



J Dae- 1 = det-i/2 f^l^ _ (^49) 



If the operator is not diagonal, the result is still the same. 

l)a{v) _ 

\/detg 



J DaC- ~ dxdva{x)g{x,y)a{y) _ iV^T (150) 
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Finally, consider 

J Dae^[°"'l, (151) 



where 



F[a,J]= J dxdya{x)g{x,y)a{y) + J dxJ{x)a{x) (152) 



is the Gaussian functional with the addition of a source term. 
For the l-dimensional case 



/OO PQO 
dxe-'^-'+b- = / dxe-«(^-*/2)=gi/4(6i6) 
-CSO — OO 

= eV4(f^f) J dxe-^^' = eV4(f^f)y^. (I53) 

Thus, 

j Dae-^l«"'l = exp Q j dxdyJ{x)g-\x - y)J{y)^ (154) 

where g~^{x — y) is the inverse operator (matrix) of g{x — y), 

/.,.-.(x-*(.-,)=*(x-„). (155) 

The source term in the Gaussian allows us to easily compute the functional integral of any moment of the Gaussian. 
Namely, 

j Daa{xi) ■ ■ ■ a(a;„)exp(- j dxdya{x)g{x,y)a{y)) 

(156) 



SJ{xi) 5J{x. 

Using the above relation, we obtain 



and 



n) J 



J=0 



Daa{xi)exp{— / dxdya{x)g{x,y)a{y)) = 0. (157) 



J Daa{xi)a{x2)exp{— J dxdya{x)g{x,y)a{y)) 



In fact, any odd moment vanishes. Here, we note that the potentially bothersome infinite number in the numerator 
is often absorbed in or canceled by a normalization. 

Fermionic Variables - Grassmann Algebra 

Since fermionic operators anticommute, we must introduce anticommuting numbers and functions, which are called 
'Grassmann variables'. Consider an eigenvalue problem for the fermion field operator, 

^l,{x)\<i>)=c|,{xm■ (159) 
For two such operators, a naive calculation shows 

V'i'02|0) = i^i4'2\4>) = 4'i<p2\4') 

= -i>2'4}l\4>) = -ll>2(pl\<l>) = -0201 10), 
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which gives (^i(/)2 = — <^2</'i- Also wc note that (fp' =Q following the Pauli exclusion principle. 
Let rj be such an anticommuting variable. This means that 

{7?,7?} = 0, (160) 

thus rf = Q. Let us define the derivative operator for an anticommuting variable in an analogous way to the bosonic 
case, 



(161) 



For any /(ry), the power series expansion of / must be of the form 

f{v) = a + br,, (162) 

since r;^ = 0. Therefore, <ff/drf = 0, so 



When a and b are anticommuting numbers, then 



We introduce integration by defining two following integrals : 

j drjT] = 1. (165) 

This definition is chosen so that / drj is translationally invariant and linear. It is not the inverse operation of d/di] ; 
in fact, integration acts like differentiation. For example, if rji and 772 are two anticommuting variables, then 

j dr}idr}2m'n2 = - j dr]idr]2'n2'ni 

= - J dT]2dr]ir]i7]2 

= -1 (166) 
Let's consider a space of anticommuting fimctions, ri{x), namely 



We carry out the functional integrals as in the bosonic case. From the definition of integration, Eq.(165), we see that 
the ^- functional d[r]] is 

m^Ilvix), (168) 

X 

so that 

J DrjS[r]] = 1. (169) 

Here, we choose the ordering of the dri{x)'s in Drj — Y[x drj{x) to be the opposite of the ordering in the functional so 
that no extra minus signs are introduced in doing the integrals. Similarly 

5[ri-cj]=^Hx)-a{x)), (170) 
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then, 

/ D7j5[rj - a]F[n] - / l[dv{x)l[{v{x) - a(a;))F[r?] 

= F[a]. (171) 

Note that F[rj] is placed on the right side of the delta function ( See p. 28, Gen. Ref. 2. ) not to spoil the ordering 

convention discussed above. 

<Problem> Prove the above equation using one component delta function S[ri — a] = r] — a. 
Now consider the Gaussian integral, 

Dr]exp{- j dxdyr}{x)g{x,y)ri{y)). (172) 
g{x,y) = —g{y,x) in order to get a non-vanishing result. For the 2- dimensional integral, we have r/i and r]2 so that 

= J dTixdme-'^"'^'^"^ (173) 

where we used r?26'2i»?i = -'n\92im = r]i9i2m- 
Expand the exponential in a power series, 

g-2r,isi.r,2 = 1 _ 2rjm2V2 + 2(r?i5i2r?2)^ + • • • 

= 1 - 2r]igi2m- (174) 

Finally, we obtain 

= j dr]idr]2{l - 2r?i5i27?2) = 2512. (175) 
Since for an antisymmetric 2x2 matrix, detg = —512521 = (512)^, we rewrite 

j dTixdme-'^'"^'^''^ = 2^/d^g. (176) 

<Problem> Show that for the 4-dimensional case, 

/.,..^.,..-<— — . 

23 

= ^(512534 - 513524 + 514523) = 4-v/det5. (177) 

We can show by induction 

j dm--- rf?7jve-2(''isi=''=+-+''"-i3''-i'"'"'' 

= (2)^/2^/^ (178) 

We assert that in the N ^ 00 limit 

j Dr/exp(- j dxdyr]{x)g{x,y)r]{y)) = (^2)°°^^^. (179) 
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Often, it is necessary to consider the fields to be complex. For bosonic functional integrals, the generalization is 
straightforward. Let a{x) be an ordinary complex function with a{x) = {ai{x) + ia2{x))/y'2 where ai{x) and 02(0;) 
are real valued. 

If we switch variables from {a*, a) to (ai, m2), the Gaussian integral can be readily performed to be 



J Da*Daexp{— J dxdya* {x)f{x,y)a{y)) 



DaiD{ia2)exjp{- j dxdy{ai{x)-ia2{x))f{x,y){ai{y)+ia2{y))) 
det/ ■ 



(180) 



In the second step, /(x, y) is assumed symmetric and can be diagonalized so that the integral gives two copies of Eq. 
(150). 

The fermionic Gaussian functional integral over complex anticommiiting functions requires a little more attention. 
Let 77 be a complex anticommuting variable and write = (■Ci + i£,2)/^f^ where ^1 and ^2 are real- valued. Then 
ri* = (a - ii2)/^/2 . 
Let 

j dr]r]=l , j difrf = 1, 

j dri = j dn* = j dm* = j dr]*r] = 0. (181) 
This means that drj = {d^i - id^2)/V^ and drj* = {d^i + id^2)/V^ ■ 

<Problem> Show that using Eq.(181), drj = (rf^i - id^2)/V2 and dr)* = {d^i + id^2)/V2 . 

Thus, we have drfdrj = —id^id^2- The 1-dimensional Gaussian integral is simply ( m : an ordinary commuting 
number) 

J dr]*driexp{~r]*mr]) = J d-q* dri{l — r]*mr]) = m. (182) 
The 2-dimensional Gaussian integral is 



/ 



dr]ldr]2dr]idr]2exp - ^ VimiVj 
= J dr]ldri2dr]idr]2 

4( » + ■ 



1 - (?7imii?7i + r]lmi2r]2 + r?2'^2i??i + ??2"^22??2) 



dr]ldri2driidri2(ri2rii'n2m{miim22 - mi2Tn2i)j 
= —detniij. (183) 
In order to avoid the annoying sign problem, we define 

Dr]*Dr] = Y[dr]*{x)drj{x). (184) 

X 

Therefore, we should order dr]ldr]idr]2df]2 • • • dr}^dr]N- Then, we have 



/ 



DrfD-qexp ^- j dxdyrf {x)g{x, y)v{y)j = detg. (185) 



The hyperbolic fermionic Gaussian integral J Drj* Drjexp (+ / dxdyrf {x)g{x,y)rj{y)) has no well-defined bosonic ana- 
logue. In order to get a well-defined limit, we should reverse the order, so that Drj* Drj = Y\xid'n{x)dri* {x)) ■ Thus, 
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/ 



Dr]* Dr]exp (^+ j dxdyr]*{x)g{x,y)r]{y)j = detg. (186) 



Functional integral techniques using the above results have emerged in last two decades as a powerful tool for study 
of many-particle systems. Unfortunately, this topic is out of the scope of this lecture. Instead, we refer to Gen. Ref. 
2 for a detailed treatment on this topic. 



B. Functional Schrodinger Picture -Schrodinger Picture in The Second Quantized Formalism 

We have learned that the Heisenberg picture method provides a convenient tool to study the many-particle system 
through Green's function. However, when the mutual interaction is strong, perturbative approaches are known less 
than reliable. Also when one studies systems such as nonequilibrium many-particle systems and evolution of a system 
in which time variation is important, Green's function method may not be a good starting point. 

In this connection, field theorists have been studying the functional Schrodinger picture field theory to provide an 
alternative theory to the Green's functions formalism. 

In quantum mechanics (first quantized system), coordinate and conjugate momentum are represented by qi{t) and 
Pi{t). In second quantized system, the field tp{x,t) and its conjugate momentum are used to represent the infinite 
degrees of freedom. 

Bosonic Theory 

We now work in a coordinate Schrodinger picture but with an infinite degree of freedom. In analogy to x\x) = x\x) 
in quantum mechanics, we consider 

^{xM=,l>{xM), (187) 

where ^(x) is the eigenvalue function corresponding to the field operator tp{x). |^) represents an eigenstate. The 
coordinate representation of the state now time dependent, is the wave functional, 'if[4']. 

^[^] = (188) 
In analogy to the quantum mechanical commutation relation [qi,Pj] = iMij, we require 

[(A(f ), 7r(x')] = iM{x - x'). (189) 

This requirement is easily satisfied by 

7r(x) = --^. (190) 

Compare this expression with the corresponding quantum mechanical one P = f ^ • 
In terms of the coordinate basis, 

(.^'|7r(f)|<^) = J^(5[<A -./.']. (191) 
With this prescription, the Schrodinger equation 

th^^\<b) = H\-^) (192) 

turns into a functional differential equation 

in-^cl>,t]=H[<l>,^—M(t>,t]. (193) 



In a time-independent case. 



H[<f>,^-^]m=Em, (194) 
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with t] = e-*-^*/^«'[(^]. 

As an example of the bosonic formahsm developed above, the free scalar theory will be studied. As it is a convention 
in the field theory, we set 7l = c = 1. The action for the free scalar theory is given (See Eq.(121).) 

= ^jd^x {d''(j>{x)d^(j){x) - mV^(a;)) • (195) 
The corresponding Hamiltonian is defined by the Legendre transform of the Lagrangian 

H = J (fx (7r(f , t) ■ dt(f>{x, t)-C). (196) 

Compare this with the quantum mechanical relation H{q,p,t) = ^^q^p^ — L{q,q,t). Here, we note that the field 
momentum conjugate to ^ is 

Using these relations, the Hamiltonian is written as 

H=^J d^x{n^ + I V<^p + m^(P^). (198) 
Using Eq.(190) for 'it{x), the time-independent Schrodinger functional equation is expressed as 

I j ^'"^ {'WW) + ^'^"^l' + m2</>2)*[<^]) = Em- (199) 
This equation can be solved exactly since there is no coupling term. 



<Problem> Solve the functional differential equation (199) using a trial functional for the ground state, \l/o[</'] = 
7/exp(— G[(^]). Show that the ground state energy is given by i^o = ^ / d^kuJk^^{Q). (See Chapt. 10, Gen. Ref. 1. ). 

Fermionic Theory 

The functional Schrodinger picture for fermionic many-particles is more complicated than the bosonic case because 
we should handle anticommuting Grassmann variables. 

In order to deal with the situation, wc examine the nature of the fermionic field theory which provides the anti- 
commutivity of the field operators. In the fermionic field theory, one can use either hermitian or non-hermitian field 
operators (For a somewhat deeper discussion, see Appendix A.). 

It should be noted that for neutral (chargclcss) Majorana fcrmions. hermitian operators should always be used, 
whereas for charged spin 1/2 particles, both hermitian and non-hermitian operators are acceptable. This comes from 
the fact that non-hermitian operators can always be written as a sum of two hermitian operators. For details of 
discussion, we refer to Appendix A and references listed there. Here, we just present the results. 

When one uses hermitian field operators, the following equal-time anticommutation relation is obtained 

{V'a(a;),i/'fc(a;')}ccO=xO' = 5ab5{x- x'). (200) 
Here, V'^ {x) = il>{x). When non-hermitian field operators are used, we obtain the equal-time anticommutation relation, 

{tpa{x),ll>l{x')}^o^^o' =5abS{x-x'), (201) 
{iJa{x),Mx')}xO=a:0' = = {ijl{x),4{x')},o^,o' (202) 

where ipHx) i^ix)- 

As mentioned above, since non-hermitian field operator can be obtained through a linear combination of hermitian 
operators, it is convenient to start with hermitian operators. Since, in this case, the anticommutation relation, 
Eq.(200), indicates that the field operator itself corresponds to the conjugate momentum (with some constant factor). 
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it is not possible to obtain a simple prescription as in the bosonic case. However, they can be represented by a linear 
combination of their eigenvalues and the functional differential operators ; 



(203) 



where 0a{x)* = Oa{x) and {9a{x), 9b{x')} = 0. It must satisfy the equal-time anticommutation relations, Eq.(200) and 
this imposes a condition for the constants a and /?, a/3 = |. 



<Problem> Show that a/3 = ^. 

Thus, a simple choice for the hermitian fermion field operators is to take \a\ 
equivalent representations 



i^a{x) = 



1 

71 



Oa{x) + 



and 



s/2i 



Oa{x) - 



%{x) 



which gives two 
(204) 

(205) 



The representation for the non-hermitian spinor field operators, which are charge eigenvalues, can be constructed 
combining two hermitian field V'l and ■^2, 



We also define a complex Grassmann variable 6{x) ; 

9{x) = ei{x) + ie2{x). 

Using Eq.(204) for both and V2, 



(206) 

(207) 



we have 



V'2 

I 



1 
1 

71 



92 + 



s 

_s_ 

592, 



(208) 



9{x) 



'{x) 



56* {x) 

5 



5e{x) 



(209) 



This prescription was first proposed by Floreanini and Jackiw [1] . 

Another way to obtain a representation for the non-hermitian field is to use Eq.(204) and (205) together with a 
real Grassmann variable 9 {9*{x) = 9{x)). 



7l(' + ^ 



From these expressions, we obtain 



■^1 



= 9{x), 



(210) 



59{x) ' 



(211) 
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which was first proposed by Duncan, Meyer-Ortmans and Roskics (DMR) [2]. It should be noted that, in this 
prescription, the hermitian conjugate of 9{x) is not given by the transpose of 6{x) but 



[^(f)]t 



(212) 



It can be shown easily that the above two representations Eq.(209) and (211) both satisfy the anticommutation 
relations Eq.(201) and (202). For charged particles including electrons, any of the above two representations can be 
used. We show, in the next chapter, examples of applications for both representations to physical problems. However, 
for neutral spin-i particles, Eq.(204) or (203) with a/3 = ^ should be used. 

Using the above prescription, we can readily write down the functional Schrodinger equation for fermionic particles. 
Especially when the Hamiltonian is time-independent. 



i^[V'^V]|*) = H 



71 (^* + ^)'7i(^^^ 



I*) 



with the Floreanini-Jackiw prescription. When the DMR prescription is used, we have 



S 



(213) 



(214) 



In the functional Schrodinger picture representation of fermionic many-particles, the physical information is given 
by the state functional {9,t\'i/) and also by the scalar product between states, 



(*l|*2) = (*2|*l)*. 



(215) 



Unlike the bosonic case, the dual to the state ket, here, is not a simple hermitian conjugate. In order to investigate 
this problem, we first consider a system of one degree of freedom described by annihilation and creation operators, 
and a, which satisfy 



{a,a^} = 1, 

{a, a} = = {a^, o^}. 

The Hilbert space of the system consists of vacuum and 1-particle states, |0) and |1) ; 

a|0) = 0,a|l) = |0) 
at|0) = |l), at|l)=0 



(216) 



(217) 



We now want to represent the Hilbert space by functional \l/[6'] of real Grassmann number 9, which satisfies 99 = 0. 
The most general functional of 6 may be written as 



^[e]=a + f30 {a,f3 = constants), 



(218) 



which implies that the Hilbert space is described by the two component basis {1,9). Thus, a natural choice for the 
^-representation is to take 



|0)^1, |1)^^, 

which imply 

(6l|0) = 1 and (6'|1) =9. 
Prom the above relation, we obtain the transformation function 

(^'|^) = (^'|0)(0|^) + (^'|l)(l|e) 



(219) 
(220) 

(221) 



This relation shows that the ^-basis is not orthonormal. With the above relation, we now define the ^-representation 
of the dual state, 
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(*il^-2) = J DeDe'{^i\e'){e'\e){e\^2) 

= j De^i[e]^2[e], (222) 

where *2[6'] = {0\^i)- Thus, the dual state *[6»] of *[6»] is defined by 

^e] = j De'{^9'){e'\e) 

= j D6'<i>^[6Y'\ (223) 

where is the Hermitian conjugate of \E'. This can be easily generalized to the case of infinite degree of freedom as 
follows 

^[9{x)]=J D0'(*|0'(x))e^'(^)^(*\ (224) 

where the variable 9 can have components in the internal space and the integral convention is used for the repeated 
variables x. 

As an example, we consider the dual state of a Gaussian functional 

vI/[6'] = ei^"^^^'^"''^^'^^''^^^ (225) 

It can be shown easily that the dual is given by another Gaussian functional (See Appendix B.) 

^[6] = C'e'^''^^''^"''*^'*^^''^^^ (226) 

where C is a constant independent of 9 and G is defined by 

G=[G^]-K (227) 

In the case of the Floreanini-Jakiw representation where complex 9 is used, the above prescription for the dual becomes 
slightly modified to become (See Appendix B.) 

§[6», 9I] = j D9' D9'^-^\9' , o''f)e^'Hs)e{S)-oHs)s'{S)^ (228) 

where 9^ denotes the transpose of the complex conjugate of 9 . However, the dual of a Gaussian functional 

5.[0,0t] =e^tGe (229) 

has the same form as in the hermitian case of Eq.(226), 

^9,9'^] = Ce^^^\ (230) 

where G = G'^ ^ . This simple property of the Gaussian functional renders the functional Schrodinger picture useful 
for computation of physical quantities. 

In the next chapter, we show some simple applications of the formalism developed so far. 

III. APPLICATIONS TO MANY-BODY PROBLEMS - VARIATIONAL APPROACH 

Many-body theoretical schemes, which arc essentially non-rclativistic field theories have been all based on the 
Heisenberg (interaction) picture, as we can see easily from the conventional Green's function method and the functional 
integral method. Only, very recently the functional Schrodinger picture method has been applied to many-particle 
systems with some successes. 

In this chapter, the functional Schrodinger picture formalism developed in Chapter II will be applied to many- 
particle systems within variational approximation. 
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A. Electron Gas and The Hubbard Model 



Electron Gas 

The Schrodinger picture formalisni provides a convenient basis for variational approximation, since, with a trial 
functional {u^ ,u\^), the ground state energy of a many-particle system can be calculated easily. 
As an example, we first study an interacting electron gas system [1]. The Hamiltonian is given by 

+ IJ2 [ d'xd'yv{x, y)i^i{S)i^l{y)My)M^)^ (231) 

a,/3 '' 

where a and /3 are spin indices. 

For field operators ip^ and ip, we use the Floreanini-Jackiw formalism, 

ip{x) = {u{x) + 



\/2 V 



x] 



*tW = ^("'W + 4=y)- (232) 

Wc note that {'ipi(x,t),%l)\(x' ,f?)^ = SijS{x — x') and u and arc anticommuting Grassmann variables. Substitution 
of field operators into the Hamiltonian turns the many-particle Schrodinger equation into a functional differential 
equation, 

H[il/V2){u + 5/Su^), {1/V2)iu^ + S/Su)]\-^o) = E\^o)- (233) 

We now calculate the ground state energy of the above many-particle system using variational approximation. A 
natural choice for the trial functional is the Gaussian, which becomes an exact solution to non-interacting many- 
particle systems. 



I^'o) = exp ^ / / d-^xd^yu'l{x)Gafj{x,y)ufj{y) 

= ex.p{u\GABUB), (234) 
(*o| =expij2 / d^xd^yul{x)Gaf3{x,y)u0{y) 

= exp(u^GABWB), (235) 



where G = (Gt)"! . Here, we used Eq.(230). 

The expectation value of the Hamiltonian is given by 



= hAB 



(^ol^o) 

«'o|i/']iV'b|*o) 



(*o|«'o) 



where 



1 (^oIV'1v':LV'bV'a|^o) , 



hAB = -^6(x~y)6af}y^ (237) 

VAB = Vix,y) = -^. (238) 
\x-y\ 
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The normalization constant (^'ol^'o) is calculated easily 

= Bets, 



where S — G + G and Eq.(186) was used. 

In order to calculate (\1/o|^]i^b|*o)) we first consider 



V'aV'bI^o) = ^ ("^ 



' SuaJ 
\- Sab 

+Gba + GscwcMr,GDA)|*o)- 



= ^{'^a'^b + Sab + ubu^jGca + u^Gbcuc 



Thus, 



(*o|V'aV'b|*o) 



Also 



(1 + C?)ba(*o|*o) + (1 + G)bc(1 - G)ba(*o|wW|*o) 
(*o|wLwc|*o) 



d 



-(*o|*o) 



d 



TtetS 



Finally, we have 



where 



dSoc dSoc 
55l,Det5 = 55l,(*o|*o). 



^BA = [{I + G)S-\I + G)]ba 

= [{I + G){G + G)-\I + G)]ba. 



Here / is the 2x2 identity matrix. 
Similarly, we obtain 

i^olijlijli^Bi^M^o) ^ 1 

(*o|*o) 4 

<Problem> Derive Eq.(244). 

Collecting the above results, we have 



[Qaa^bb — ^ba^ab]- 



2m 



S{x-y)V^Q.aa{x,y) 



a,l3 •' 

- d^xd^yV{x,y)n/3a{y,x)nai3{x,y). 
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Here, we note that the second term corresponds to the direct term and the third term to the exchange term, respec- 
tively. 

When the system is homogeneous, it is convenient to use Fourier transformations of 0^/3 and V{x,y) . 

-±^e^'^<'^-y)u^0{k), (246) 
^e^H^-y^k). (247) 

V{k) = 47re^/|fcp for the Coulomb gas. 

These prescriptions allow one to express the expectation value of the Hamiltonian in fc-space. 

{H) = ^VgY^ j dKeo{k)wUk) 

a 

-IVgY.! dKdK'V{k - k')uf3„{kPaf3{k'), (248) 

where eo(A:) = dK = and Vg is the volume of the system. 

In order to simplify the above expression, we introduce the particle number N, 



N : 



\^9Y. j dKu^aik), (249) 



which follows directly from Eq.(242). 
Thus, 



with 



= ^^Tr[H + IvimVVg, (250) 

haoik) = eo{k)6a0 -\ j dK'V{k - k'^^fiik'). (251) 

In order to obtain the ground state energy, it is necessary to take a variation on {H) with respect to g or g, which 
is the Fourier transformations of G or G respectively. 

^=2^./^..^- 0. (252) 

Using u) = {I + g){g + g)~^{I + g), we obtain 

^(/-,)M/ + ,) = 0. 

Similarly, we can obtain the condition for g. Altogether, 

{I-g)h{I + g)=0, {I + g)h{I-g) = 0. (253) 

From Eq.(253), we obtain 3 types of solutions, namely, g = — /, g = I and the last type g / ±7. 

g = —I is a trivial solution which gives lu = 0. Next wc consider g = I. It is known that it gives a trivial solution 
in the free Dirac theory and the Gross- Neveu model [2, Gen. Ref. 5]. However, it will be shown that g = I gives 
a Hartree-Fock solution in the electron gas model. We believe that this interesting phenomenon arises due to the 
difference in the nature of the ground states in the electron gas from that in the Dirac theory. In the Dirac theory. 
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the vacuum is the Dirac sea with an infinite energy which should be renormahzcd away, whereas in the electron gas, 
the vacuum ground state is the fermi sphere which is a very much relevant quantity in condensed matter physics. 
When g = I, we obtain u = 21. Therefore, at T=OK, the particle density Uapik) is given by 

uJo.f){k) = 2Sa.fje{kF - fc), (254) 

where kp is the fermi wavevector and 9{kF — k) = I for k < kp and for > kp. 

When this result is substituted in Eq.(250), we obtain the variational ground state energy in the Gaussian functional 
as follows, 

(H) = {2S + l)Vg [ dK[eoik) + S(fc)]0(fcf - k), (255) 



where S represents the spin. The self-energy T,{k) is defined as 



E(fc) = lnV{0) - I I dK'V{k-k')9{kF - k). (256) 



Here n = N/Vg is the electron density. We note the above result is the standard Hartree-Fock result (Eq. (10.22) in 
Gen. Ref. 3). 

Hubbard Model 

Most of the theories of strongly correlated electron systems begin with the Hubbard model [3]. Therefore, it is a 
natural choice to be studied by the new many-body theory via the functional Schrodinger picture. 
The Hubbard model Hamiltonian is given by 

H = - '^eijcl^Cja + h.c. + U y^n^in^ 

ija i 

= -Yl ^abc\cb + h.c. + ^UY^ TABnAUB- (257) 

A,B AB 

Here, the first term represents the hopping term, the second the hermitian conjugate of the first term, and the third 
the on-site Coulomb repulsion. We used the notation A = (i,a), B = {j,l3), Iab = (ij^ais, and Tab = — Sa^)- 
i and j represent the lattice site and a and /3 the spin indices. 

Since, now, we are dealing with discrete lattice sites, we should generalize the Floreanini-Jackiw formalism to the 
discrete case by writing 

Here, ua and u\ arc Grassmann variables which satisfy {ua,u^b} ~ ^ab ■ The Gaussian trial wavefunction assumes 
the same form as Eqs.(234) and (235), 

(^ol = e^'^^GABUB^ (259) 

but with the summations on i and j instead of integrals on x and y. Following the same steps, Eqs.(236) (244), we 
obtain 



A,B 



+ ^UY^TABi^AA^BB - ^bA^Ab), (260) 



A,B 



where n = {I + G){G + G)-\I + G) . 

In terms oi i,j and a, /3, the expectation value becomes 
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ijal3 

+ X [(^^ji)c<a(^M)/3/3 ~ (^M)/3a(^M)c</3] 

= — ^ ^ Cij {^ij)aa -\- h.C. 
ija 

+ gC^ ^ [(fiji)aa(f^M)/3/3 - (f^M)/3a(f^M)a/3]- (261) 

Here, we assume the translational symmetry and consider only the nearest neighbor hopping to. Then, the fc-space 
expression is given by 

{H) = -to^^k[(^k]ac, 

+ ^ [{^k)aa{i^k')l3(3 - {l^k)l3aii^k')al3], (262) 

kk'aP 

where. A'' is the number of the lattice sites, u) = {I + g){g + g)~^{I + g) with g the Fourier transformation of G. ^k is 
defined 



^^=Y^^k.r (263) 



where r is the vector to the nearest neighbor. 

In order to study various magnetic ground state, we rewrite {H) using 



k.OL 



2N 

k 



which turns out to be the particle number per site. {H) is now rewritten as, 

{H) = \Y,Tv[h^]-\uNn\ (265) 



where 



hk = [-2to7fc + ^^]^ - J ^ E ^fc' • (266) 

k' 

Taking a variation on {H) with respect to g as done in Eq.(252), we obtain 

{I-g)h{I + g) = Q, (267) 

which has the identical form for the electron gas model. 

As we did for the electron gas model, we first consider g = I. In such a case, as in Eq.(254), we readily find that 

u = 29{kF - k)I. (268) 

Using this expression and the relation ^^.^ 9{kF — k) = Nn, we obtain 

3 

hk = [-2io7fe + -^Un]I, 

and 
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hu = 2[-2io7fc + ^Un]e{kF - k)I. (269) 

Finally, {H) is given by 

{H) = -2to ^kO{kF -k) + ^n^NU, (270) 

ka 

which is the well known paramagnetic mean field result of the Hubbard model [3,4] 
In order to study magnetic ground states, we introduce the following order parameter 

fe 

which is the difference between the up-spin and the down-spin electron divided by the number of the sites. 
The expectation value of the Hamiltonian can now be rewritten in terms of n and m, 



{H) = \Y. T-^M - luNln" - m\ (272) 

where 



2 

fe 



h = [-2<o7fe + \un\ao - ^Umas- (273) 

<Problem> Show that Eq.(272) is identical to Eq.(265). 

Again taking a variation with respect to g, one obtains 

{l-g)hil + g)=0. (274) 

The above equation can be solved by expanding g and h in terms of the Pauli spin matrices. However, since wc seek 
a solution relevant to ferromagnetic ground state, we assume that g and h have ctq and (73 components only. Thus, 
[g, h] = and (/ — g^)h = 0. Therefore we have 

g = ±a3 or ±7. (275) 

Since g = I has been studied already, g = ±173 will be studied here. 

9 = <y<J3, (276) 

a=i \ up- spin _ 

I —1 down — spin ^ ' 

Now, in general w = (Tq + eras and can be written 

Lo = 9{kFa - k){ao + (TfJs) (278) 
at zero temperature, hco in Eq.(272) is now expressed as a diagonal 2x2 matrix 

+ <^)iho + WkF, - k) \ 

V - (^)(ho - h3)9[kFi - k) J ' ^ ' 

where ho = —tajk + \Un and /13 = — ^Um . This expression allows to write the 

{H) = Ek^e{kF^ -k) + Y EkiOikFi - k) 

k k 

-^C/iV(n2_TO2), (280) 

where 
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Ek-[ = -2io7fc + ^U{n - m), 

Eki=-2mk + \u{n + m). (281) 

The above result clearly shows that the splitting due to the mean-field magnetization. 
It can also be written in a more familiar form using ^^.^ \Un = ^UNtt?. 

{H) =Y^Eka + \uN{n'' + m^), (282) 

where 



ka 



Eka = -2to7fe - ^(^Um, (283) 

which is the ferromagnetic Hartree-Fock result [4]. 

A more interesting case is provided by the antifcrromagnctic ordering. For simplicity, only one dimensional case 
will be considered. The order parameter in this case is the staggered magnetization, 

mi = (-l)'m = me'^'', (284) 

where Q = ir/a and r = la. This is readily reduced to the paramagnetic state by setting m, = 0, and to the 
ferromagnetic state by setting Q = 0. The Hamiltonian can be expressed in terms of this order parameter as follows, 

{H) = -\n j dKTi{hkUJk) + \um^ + ^Un'', (285) 

where 

hk = €kao + Mgas. (286) 

Here, dK = dk/2'K and Mq = ^Ume'^^ 

For the antiferromagnetic ordering, the Brillouine zone should be folded in half. This zone folding can be carried out 
by separating the fc-contributions into two parts 

{H)=-In [ dxf ^Q''^]( )+}.Um' + }.Un^ (287) 

where the prime on the integral signifies the zone-folding. This Hamiltonian is diagonalized under the nesting condition 
€k+Q = —Ck- Now, we take a variation in g and obtain (ctq — g)h{ao +g) =0. Among the solutions, we choose g = ao, 
because the antiferromagnetic gap occurs due to the nesting and the zone-folding and, thus, requires no further 
symmetry breaking in terms of spin variables. Finally, we obtain 



(288) 



which gives two antiferromagnetic bands degenerate with spin indices [4]. The paramagnetic ground state is readily 
obtained by requiring m — 0, and the ferromagnetic state is obtained by setting Q = in Eq.(286). 

When the on-site Coulomb interaction U is large (U ^ t), doubly occupied sites are energetically very expensive. 
Thus, the restricted Hilbert space now consists of configurations made of empty sites (holes), and up and down 
spins. The effective Hamiltonian which describes this situation is called t-J model and has been studied intensively in 
possible connection with high Tc superconductivity and 2-dimensional antiferromagnetism [5]. Functional Schrodinger 
picture approaches coupled with the slave-boson and the slave-fermion techniques have been shown to produce correct 
mean- field results on the 2-dimensional t-J systems [6]. 

The above results shows versatility of the functional Schrodinger picture in obtaining various ground state energies 
within the Gaussian approximation. However, at the same time, one may notice that the variational method coupled 
with a Gaussian trial functional invariably leads into mean-field or Hartrec-Fock results. This has been early noted 
that in the field theory and the mean-field results are often called Gaussian. In the later sections and chapters, efforts 
going beyond the Gaussian approximation will be described. 
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B. BCS Superconductivity 



The phenomenon of superconductivity is a rare example of a simple field theoretic many-body theory being able 
to explain most of the essential features. Since Bardeen, Cooper, and Schricffer(BCS) introduced their variational 
theory in 1957, several alternative methods including the canonical transformation technique and the Green's function 
method have been successfully applied to superconductivity [Gen. Ref. 3] . Therefore, it appears imperative that any 
new many-body theoretical technique should stand the test of the BCS superconductivity. 

In this section, the functional Schrodinger picture formalism is applied to the BCS superconductivity [1]. Here, 
instead of the Floreanini-Jackiw formalism, we employ the DMR formalism, Eq.(211) to solve the BCS equation. It 
was shown that the DMR formalism require a non- conventional definition of the hermitian conjugate. However, it 
will be shown that in order to obtain the gap equation, we do not require the dual of the ground state. Therefore, 
no complications due to non-standard hermitian conjugate arise. The same problem can also be solved using the 
Floreanini-Jackiw formalism [2]. But it will be left as a problem. 

The model grand canonical Hamiltonian for the BCS superconductivity is given by [Gen. Ref. 3] 

H^J2j d'^^ii^) (-^V^ - Mx) 

+ J^E / c^'a:^i(2^)»^SV'M^)#(^)*^?iV'a(:?), (289) 



where fi is the chemical potential and 



-i 

1 



is introduced to insure the s-type pairing. In the process of factorizing the quadratic interaction, the usual Hartee-Fock 
terms are neglected. Eq.(289) is, then, reduced to 

a •' 

al3 '' 

+ V^-^, (290) 



when^(f) = -{ti^/2m)V^5{x'-x) — ji. x' denotes a coordinate infinitesimally larger then x, which ensures the proper 
operation of ^{x) on tjja{x). The gap function is defined by 

A = -g{i>a{x)ia'^^^7p0{x)) 
= -5'(V'T(^)V'i(^)) 

= ff(Vi(f)VT(^))- (291) 
In the DMR formalism, the fermion field operators are expressed as 

where ip{x) is a real Grassmann function. With the above prescription, the time-independent Schodinger equation is 
given by 
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(293) 



where $[1^3] is the wavefunctional which satisfies the functional differential equation. 

Here, it should be noted that we do not need the troublesome dual of since it is not necessary to evaluate {H). 
Power counting indicates that G[if] should minimally be quadratic in ip so that the second functional derivative of G 
yields a number while the square of the first derivative is still quadratic in 

Therefore, we write G[ip] as a quadratic form in ^p ( For the power counting technique, see Chapter. 10, Gen. Ref. 
1.). 

In the mean field theory of the electron gas and the Hubbard model in the previous section, G has a form, V 
However, the DMR formalism employed in this section requires G to be written in the form ipFip, which reflects the 
pairing. With this prescription, G[(p] is written 



a/3 

We note that Fap{x, y) must be symmetric in coordinate space and antisymmetric in spin space so that ^[f] represents 
the s-wave pairing. Substituting Eq.(294) into Eq.(293) and equating terms with same order, we obtain two equations. 




(294) 




(295) 




(296) 



Eq.(296) can be written in a simple matrix form 




where K = Fa^^\ It is convenient to solve Eq.(297) in momentum space. For this purpose, we define K{k), the 
Fourier transformation of K{x, y), 




i=0 



(298) 



where dK = d^k/{2'K)^. This Fourier transformation reduces Eq.(297) to a set of coupled equations. 



3 




2^kKi + AiAKoKi = 



(299) 



where — Ek — n = h'^k'^ /2m — ^. 

The above equations are readily solved to yield 




(300) 



Now F{x, y) is given as 




(301) 
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where Ek = ^/^^-\-\A^ represents the excitation energy. The wavefunctional which satisfies Eq.(293) is given by 

*M =expQ I dK ^!^^[^^{k)^i{-k) - ipi{k)ip^{-k)\^ . (302) 

It is obvious that the above wavefunctional has the correct BCS state features. Also the ground state energy is 
obtained from Eq.(295) and given by 

E = V J dK i^k- Ek) + (303) 

Here, A is assumed real. Since the characteristic energy range in forming the Cooper pairs is hwc, the ground state 
energy per volume is 

cs = E,/V = N{0) / deU- ^/eTA^) + —, (304) 
Jkw^ ^ ^9 

where hwc is the phonon energy and -/V(O) the density of states at the Fermi energy . 
The gap equation can be obtained by minimizing Eq.(303) with respect to A, 

A=f/7?i^^, (305) 
2 7 (27r)3^^2 + A2' ^ ' 

which can be reduced to the standard BCS result [3] 

^ = sinh[l/Ar(0)5]- ^^0^) 

It is shown that the functional Schrodingcr picture method via the DMR formalism can be used to construct a consis- 
tent theory of the BCS superconductivity. It should be noted that we could avoid the problem of the nonconvcntional 
hermitian conjugate in the DMR representation since one does not have to calculate the expectation value of the 
Hamiltonian. It is because that the factorized form of the BCS Hamiltonian is only quadratic and, thus, can be solved 
exactly without resorting to the variational process. 

Since, the Floreanini-Jackiw formalism is also a valid representation for fermions, it should be possible to solve the 
BCS superconductivity problem using this formalism. This is done in the literature and will be left as a problem. 

<Problem> Obtain the gap equation using the Floreanini-Jackiw formalism (Reference. 2). 



C. Dilute Bose Gas and Bose-Einstein Condensation 



It has been shown that, in Section II. B, the functional Schrodinger picture formalism assumes an easy and 
straightforward form in the bosonic scalar field theory. However, when this formalism is applied to boson particles, 
occurrence of the Bose-Einstein condensation makes the calculation less straightforward. 

In this section, the bosonic functional Schodinger picture is generalized to include the Bose-Einstein condensation. 
Also, we try a shifted Gaussian trial functional in order to reach beyond the Gaussian approximation [1]. 

The standard model Hamiltonian for an interacting Bose gas is given by [Chap. 10, Gen. Ref. 3] 

h^k"^ U 

-'^ X] ^;^"fc"'= + ^ 4i42«fe3«fe4^fei-Ffc2,fc3-Ffe4- (307) 

k kik2k3ki 

The constant matrix element U can be determined by requiring that H correctly reproduce the two-body scattering 
properties in vacuum. We assume that the scattering length is much less than the inter-particle spacing n^/^, where 
n = N/V. Therefore, the validity of the present calculation is restricted by the condition na^ <C 1. Under this 
assumption, the zero momentum operators ao and a\ become nearly classical due to occurrence of condensate. The 
interacting part of the Hamiltonian is rewritten explicitly in terms of ao and Oq as 
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Hint = ^alalaoao + |^2(o|.ajo/sao + aL^aJa_feao) 

+ ala'_^aoao + alalaka-k^ 



(308) 



where only terms of order of A^q , A'o and ^/Nq have been retained. The three-zero-momcntum operator terms do not 
exist, because momentum conservation makes the fourth term also have zero momentum. It should be noted that 
the last term, which is generally neglected in text book calculations originates from the interaction of particles with 
condensate, thus representing the condensation fluctuation. 

In order to handle the contribution from the last term, we introduce a new variable jk, 



7fc = X]4+9S- 



(309) 



We approximate as a c- number. Replacing the operators uq and Oq by \/No, Hint becomes 



Hint — 



U 
2V 



U 

2V^ 



u 



k^O 



The number operator satisfies the relation 



(310) 



(311) 



k^O 



Using this relation and keeping up to the V A'' terms, the model Hamiltonain becomes 

H = ^n'^VU + i ^ [{el + n[/)(4a, + al.a.k) 



kjtO 



+ n[/(4aLfe + aka-k) + ^ '^'^^^'^ + ^-k)^ 



(312) 



where e° = h^k^ /2m,. 

When there is no Bose-Einstein condensation, the bosonic operators can be simply written = (p{k) and aj. = 
—5/5(j){k). However, due to Bose-Einstein condensation, such simple prescription does not work. With reasons which 
will become transparent below, we introduce Bose-Einstein weighting factor, ak{k ^ 0) and express the operators. 



dk = cxkcpik), 4 = -ak 



S 



6(j){k) 



(313) 



When these expressions are introduced to the commutation relation [afe,a^] = Skk', we obtain a constraint relation 

for ak, 

namely 



akak'S{k - k') = 5, 



kk'- 



(314) 



Note that q:^(5(0) = 1 except k = 0. The Hamiltonian is now given by 



Him,- 



5(t>{k) 



\Y.'^l[-{^l + nU) 



Scpik) 



m 



S(j){-k) 
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+ nU 



5(j){k) 5(j){-k) 

m 



nil \ ^ 



s 



5(j){-k) 



Here, we assume that the ground state functional is real and positive and has no nodes. Thus we write 



(315) 



(316) 



where C is the normahzation constant. 
The Schodinger equation becomes 



fe^O 



nU 



kjtO 



SG 6G 



nU \ ^ 



6(l){k)S(p{-k) S(j){k) 6(t){-k) 
SG 



m 



Scpi-k) 



Eg. 



(317) 



The power counting technique indicates that G[(/)] should be at least quadratic in (j){k). Therefore, a simple form of 
G[(/)] can be expressed as a shifted Gaussian in ^, 



G[ct>]= j d^k[g{k)m^{-k) + f{k)4>{k% 



(318) 



where (j){k) = —4>{—k) to ensure convergence and ^(A:) and f{k) are function parameters that will be determined below. 
It should be noted that the linear term in the above expression has no contribution if the fluctuation contribution is 
not included in the calculation. Substituting Eq.(318) into Eq.(317), we obtain the following relation. 



Eg = \n'VU + ^ al5mnUg{k) - (e° + nU)] 



nil \ ^ 



k^tO 



N 



k^O 



nlJjk 



-ie° + nU)f{k) + 2nUg{k)m + i|(l - 2g{k)) 



i ^ [ - 4(e° + nU)g{k) + AnUg{kf + nU] m<t>{-k) 



k^O 



Power counting of the above equation readily yields, 

AnUg{kf - 4(e5! + nU)g{k) + nC/ = 0, 



Eg = \n^VU + Y.[nUg{k) - {e^ + nU)] 



-{el + nU)f{k) + 2nUg{k)f{k) + ^(1 - 2g{k)) = 0, 



nil \ ^ 



k^O 



k^O 

aimn-k) - 



N 



(319) 



(320) 
(321) 



(322) 



In deriving the above relations, we used Eq.(314). Function parameters are now determined to be 
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_ el + nU + E{k) 



Ik 



E{k) 



ak^/N E{k) 



where 



f{k) 



E{k) = J{el+nUY-{nUy. 



(323) 

(324) 

(325) 



Finally, the ground state energy is given by 



Eg = Eq + El + E2 

= ln'VU-lYlii4 + nU)-E{k)] 



Ikl-k 



(326) 



The first two terms, £^0 and £^1, are the standard results which can be obtained by a canonical transformation [p. 317, 
Gen. Rcf. 3]. The last term, E2, is the fluctuation contribution from the condensation. Following the standard step 
relating the scattering length to the interaction [p. 314, Gen. Ref. 3] 



m 



rr 1 



we obtain the ground state energy density 



E/V = 



2Trh an 



m 



1 



128 
15v^ 



3U/2 



(327) 



(328) 



which is the standard textbook result. Note that we have not yet evaluated the fluctuation contribution from the 
condensation. 

Before we evaluate the new correction term to the ground state energy, it is interesting to check whether the present 
formalism predicts the particle depletion from the zero momentum correctly. 
The depletion is defined by 



N-N, 



N 



^ V^/ t \ 



(329) 



fc#0 



This value is equal to 70 /-/V and can be calculated using Eq.(313), (314) and (316). We obtain for nu = (^Ofc) 

5G 



Oil 



^(0) + </.(fc)- 



5(j){k) 

l + 2alg{k)m' - 4fik)m), 



where we used 4'{k) = — (/>(— fc). 

The expectation value of <p{k)'^ and 0(fc) can be evaluated by carrying out the Gaussian integrals 



(330) 



($o|0(fc)2|$o) 



1 



and 



(<E.o|$o) Mk) 4<7(A;)2' 
($o|0(A)|$o) f{k) 



(*o|*o) 



2g{k)- 



(331) 



(332) 



<Problem> Prove Eqs.(331) and (332). 
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Using these expressions, we obtain 



nk 



-1 H 

^ 2 ' 



(333) 



where k . 

In order to obtain ak, which is introduced as a measure in the k-space to conserve the particle numbers, we use a 
simple physical argument. We note that rik should satisfy the following conditions 



lim rifc = oo and lim Uk 

k — *0 k — >C30 



0. 



Also, we obtain from Eq.(323) the following properties 

1 



5(0) 



and £(0)~0 for fc ^ 



and 



nUg{k) - E{k) for k 



oo. 



The simplest function to satisfy the above properties for Uk is given by 

nUg{k) 



at 



E{k) 



(334) 



(335) 



(336) 



(337) 



Indeed, we note that Uk has the desired properties of a very sharp peak at A; = 
The particle depletion is now calculated using Eqs.(333) and (337) to yield 



and almost zero elsewhere. 



N-No 
N 



-y 

fc#0 



nUg{k) 



k^to 



m 

'el+nU 



(338) 



1/2 



which gives the standard textbook result, | (^^^j after carrying out the k summation [p. 317, Gen. Rcf. 3]. So 

far, we have shown that the functional Schrodinger picture can handle dilute bose gas successfully even when the 
Bose-Einstein condensation occurs. Next, it will be shown that the functional Schrodinger picture scheme coupled 
with a shifted Gaussian function enables us to calculate the fluctuation contribution from the condensation. 
The ground state energy of dilute bose gas is generally given (See for example Eq.(22.21), Gen. Ref. 3.) 



E 
V 



m 



128 



na 

TT 



1/2 



+8{-n-Vs] {naf In(na^) + 0{na^) + ■ 



(339) 



We note that the first two terms are the Gaussian results discussed above. The next order correction term originates 
from the second term of Hjni, Eq.(308). However, the coefficient of the last term has never been determined. 

Below, we show that the fluctuation contribution, the last term in Eq.(326) produces a term in na'^. Indeed, 
a rough estimate of this contribution showed that this contribution would be of the order na^ [p. 318, Gen. Ref. 
3]. The fluctuation contribution, £2, carries the same ultraviolet divergence as in Ei. The divergence in Ei was 
handled through the expression of the scattering length a to order of as given in Eq.(327). However, for the 
calculation of E2, such a simple cancellation scheme does not exist, because it is now necessary to expand a up to 
U^. Therefore, we resort to an alternative cutoff procedure called minimal subtraction [2]. In minimal subtraction, 
ultraviolet divergences are removed as part of the regularization power of k from the momentum space integrated. 
Following this prescription, we obtain E2 in a dimensionless form 



E2 



2U 

'nv 



fe#0 



1 



el + 2nU e 



(340) 
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where Ec — 2nti^an^V/m . Now, it is necessary to evaluate 7fc7_fc- It is known that a dominant contribution arises 
when a particle interacts with itself, which happens when f = k + q. Therefore, 



qq' 



E 



6° + nf/ - E{q) 



-\ 2 



2E{q) 

This approximation gives a simple result for the condensation fluctuation correction E, 



E2/EC 



16 



- na 



(341) 



(342) 



Here, we note that the above result is obtained through a rather drastic approximation of Eq.(341). Also, there 
are indications that there might be na^ order terms from other sources [p. 319, Gen. Ref. 3]. A calculation which 
produces a na'^ order term from the particle-particle interaction has been reported [3] . 

In this section, we generalized the bosonic functional Schrodinger picture to the case when Bose-Einstein con- 
densation exists. Also, a shifted Gaussian functional is employed to calculate the fluctuation contribution from the 
condensation. 



D. Finite Temperature Formalism of The Electron Gas 



So far, we have presented only zero-temperature formalism. In order that the functional Schrodinger picture 
many-body theory be truly useful, it requires extension to finite temperature. In this section, a finite temperature 
formulation of many-body theory based on the functional Schrodinger picture is presented [1]. It is shown that a 
Gaussian approximation of the theory produces the finite temperature Hartree-Fock results both for the para- and 
ferromagnetic phases in a simple and convenient fashion. It is also shown that the self-consistent equation for the 
ferromagnetic splitting yields the Stoner condition for the ferromagnetism. 

A finite temperature quantum field theory using the functional Schrodinger picture has been recently formulated 
by generalizing the Gaussian approximation of the variational calculation in quantum mechanics [2]. The method is 
based on density matrix formalism in statistical mechanics [3] . 

Here, we briefly introduce the concept of the density matrix. When Icpi) is an eigenket, Ei the corresponding 
eigenvalue of the Hamiltonian H of the system and Q the partionfunction, the probability that the system is in the 
state \4>i) is {1/Q)e~^^\ Thus, the density matrix is defined to be 

p = Y^n\<Pn){(l}n\, (343) 

n 

where a;„ = (1/Q)e-^^" and P = l/ksT . 

We can readily show that {A) = TvpA and Trp = 1. Because H\<j)n) = En\(j)n), we can write 

1 ^ p-l^" 

P=QY^~^''''\^n){<i>n\ = (344) 
^ n 

where 

e-P^ = Q = Y ^'"^^ = T^e-^", (345) 

n 

and, thus. 
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(346) 



Now we regard the density matrix as a function of f3. The unnormaUzed p is defined by 

PuiP) = e-^". (347) 
In place of pu, we will hereafter write p. p{(3) obeys the differential equation 

-%=Hp. (348) 

<Problem> Using the energy representation, H\i) = Ei\i), prove Eq.(348). 

Solving this differential equation, one can obtain p(/3), F and other physical properties of the model system repre- 
sented by H [3]. In order to develop a finite temperature many-body theory, we first study non-interacting fermion 
system which has a grand Hamiltonian given by 



K = H-fiN 





2m 



iPccix) - fiN, (349) 



where N represents the number operator and p, the chemical potential. Using the Floreanini-Jackiw prescription, the 
grand Hamiltonian is expressed as 

^ = iE«-(«U5^)(". + 4). (350) 

where £,ab = -{h"^ /2m)[6{x - 27)V^](5„/3 - p5ab and A = {x,a), B = {y,l3). 
We choose a Gaussian trial density matrix 

p = exp{u\FABUB + JaUa + u\La + C). (351) 

to satisfy Eq.(348). 

Substituting Eq.(351) into Eq.(348), we obtain the following relation. 

_9J _ ^dL_dC_ 

+ {^ + F^-J^L)\. (352) 
Since the above equation is valid for any Grassmann variables u and u\ we obtain four differential equations, 

(353) 



dC 

'dp 



Ti[{I + F - JL)^]. (356) 



These equations are readily solved to yield 
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F = -coth( -/3C 



sinh ( 



-ve 2 



sinh ( 



C = Tr 



sinh ( 



— t;^vcoth f 



(357) 
(358) 
(359) 
(360) 



In the above Grassmann variables v and are introduced as constants of integrations in order to make p as a Gaussian 

type density matrix. 

Collecting these expressions, we obtain 



{PnorU = -ie^[-i^«]Det 



sinh 



v—v ' u 



^_(„t+^t)[coth(l/3|)](u+t;)^ 

where the normalization constant A'' is to be obtained by the relation, Tr[p„or] = 1- 
<Problem> Prove Eqs.(357) ~ (361). 

We note that N = Tr[p]. Carrying out the integral, we obtain 

UU 



(361) 



= eT^[-^'^«]Det 
= e^[-i^«]Det 



sinh ( -/3^ 



4 cosh ( 



where we used Eq.(185). The normalized density corresponding to the above density matrix is given by 



Pnor — 



Det [4cosh(i/3^)] 
The Helmholtz free energy is defined by 

l3F={\np^or)+P{H)- 
Substituting Eqs.(349) and (363) into Eq.(364), we obtain 

F = Nii- keTTr [in (l + e"''*)] - Tr [In 2] 
= NiJ.- keTY^ [In (l + e"^**^)] - Arin2, 

where ^fe = ^ M ^^^^d is the spin index. For the second expression, we carried out the Fourier transform. 
The total energy is calculated through 

E=(H) 



(362) 

(363) 

(364) 

(365) 



= ^ / Du'DuhAB'>pA'>pBPu 
AB •' 



Puv 



AB 
1 



-Tr <^ h 



1 - tanh ( -P^ 



(366) 
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where Hab = £,ab + M and rife = 1/ {l + e^'^^^"'*)^ . 
<Problem> Prove Eqs.(365) and (366). 

We observe the well-known results of non-interacting fermionic particles, Eqs.(365) and (366) clearly confirm that 
a finite temperature many-particle theory can be successfully formulated using the density matrix method based on 
the functional Schrodinger picture approach. 

Next, we investigate the interacting electron gas system of which the grand Hamiltonian is given by 

K = Y1 Ub^Wb + VAB^Wsi^Bi^A, (367) 

AB AB 

where Vab is any appropriate particle-particle interaction potential. The density matrix for an interacting system is 
not a Gaussian and generally quite complicated. As a first approximation, we choose a trial Gaussian density with 
a variable function, which will be chosen from minimization of the thermodynamic potential. Thus, the normalized 
trial density matrix is given by 

r 1 1 

[Pnor\i 



" Det[4coth(/jg/2)] 

^^-(u^W)[coth{0Q/2)]{u+v)^ ^3gg^ 

where Q is the adjustable parameter matrix. The corresponding normalized density operator is written in the form 



pTr[/3Q/2] , 

o - ^ f.-0i'^Q^ C369^ 

Det[4cosh(/3g/2)]^ ' ^^^^^ 



in analogy to the free particle case. 

The thermodynamic potential is given by 



/JO = (i{F - nN) 

= {lnpnor)+P{H)-0f,p{N) 

= {\npnor)+l3{K). (370) 



Substituting Eq.(369) into Eq.(370), one can obtain the following expression 



/JO = Tr[/3Q/2] - (lnDet[4cosh(/JQ/2)]) 

-/3(V^QV)+W- (371) 



The third term is calculated using the same steps used in Eq.(366). 



AB 



= I:QabJ Du^Du\ {u\ + ^) [uB + ^ j p. 

= IY^QabGba, (372) 



'ab 



where G = ao — tanh(/3(5/2). 
Similarly, we obtain 



{^^^i>) = lY.^ABGBA, (373) 

AB 

(tp^BVABi^B^A) = 2 H Vab[GaaGbb - GbaGab]- (374) 



^AB 



Collecting terms, the thermodynamic potential is expressed as 
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In 



4 cosh ( ^I3Q 



A A 

--^P'^QabGba + ■^P^S.abGba 

AB AB 

+\i3Y.Vab[GaaGbb-GbaGab]. (375) 



AB 

For a homogenous system, it is more convenient to carry out the calculation in fc-space. Thus, 



pfl = -V j dii'tr[ln[4coshy4(fc)]] 
+V J dii'tr[A(fc)tanhyl(fc)] 

+ ^VP I dKtrUk)G{k)] 



+^VI3V{{)) j dKtr[G{k)] j dK'tv[G{k')] 

~V(5 j dKdK'V{k-k')ti[G{k)G{k')], (376) 

where A{k) = (3Q{k)/2, G{k) = ctq — tanh(/3Q(fc)/2) and tr is the trace for the spin indices only. 
Introducing a constraint parameter N, 

N = ^V j dKtv[G{k)], (377) 
we rewrite the thermodynamic potential as follows 

I3fl = a N-^V j dKti[Gik)] 
-V j (iii'tr[ln[4coshy4(fc)]] 
+V J dKtr[A{k)ta.nhA{k)] 
+ ^VP J dKtr[^{k)G{k)] 

-^VP j dKdK'V{k - k')tY[G{k)G{k')], (378) 
where a is the Lagrange's undetermined multiplier. Taking variations on /30 with respect to A'' and A{k), we obtain 



a = -ViO)f3N/V, (379) 
1 

2^ 



f3Q{k) = iPiik) - a)ao - I dK'V{k - k')G{k'). (380) 



<Problem> Prove Eq.(380). 



Eq.(380) is a self-consistent equation which determines Q{k). Since all functions in the above equation are 2x2 
matrices, Q{k) and G{k) are also expressed as linear combinations of Pauli matrices. First, we consider the symmetry 
conserving paramagnetic case. In this case, the eigenvalues of G{k) are same for both spins, and, thus, Q{k) and G(k) 
must be multiples of the 2x2 identity matrix same as in the zero-temperature case. 
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Q{k) = q{k)ao, 

G{k)=g{k)ao. (381) 



By defining q{k) = Ck — l-i, g{k) can be written as 

"1 



g{k) = 1 — tanh 
2 



1 _(_ Q0(f-k-IJ.) 



2nk. (382) 



Here, Uk is the Fermi distribution function of the interacting system. The energy spectrum of the interacting system 
is obtained from Eq.(380) 



ek = el + {2S + l)V{0)n - j dK'V{k - k')nk' 



(383) 



where n = N/V is the number density. The above equation is the famiUar self-consistent Hartree-Fock equation, 
which determines the energy spectrum of the interacting electron gas system. The total energy is obtained using 
Eqs.(373), (374), (381) and (382). 

E={H) 

= \y j dKti[h{k)G{k)] 

+ ^VV{0) j dKtr[G{k)] j dK'tv[G{k')] 

-\v j dKdK'V{k-k')iv[Gik)G{k')] 
= {2S+l)V j dKeluk 

+ ^VV{0){2S +l f j dK2nk j dK'2nk' 

-^V{2S+1) j dKdK'V{k-k')Anknk' 
= {2S + l)V j dK{el + hlk)nk, (384) 

where the self-energy Efc is given by 

Y^k = (25 + l)V{Q)n - j dK'V{k - k')nk- ■ (385) 

The thermodynamic potential is obtained using Eqs.(376), and (381) ~ (384), 

0= -NkBTln2 

- kBT{2S + 1)V j dK In l + e-(^'=-'')/'=^'^ . (386) 

<Problem> Prove Eq.(386). 

As we discussed in Section III. A, there exist solutions other than the simple paramagnetic solution. In order to 
consider a magnetic symmetry breaking solution, we express Q{k), 

Q{k) = [ck - n)(To + ikCTi. (387) 
Assuming that 7fe is much smaller than Ck and expanding to the first order of 7^, we obtain 
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tanh 



-/3 [(efc - /u)cro + jkos] 



(To tanh 



1 2 
+ f73 2/37fcsech 



-I3{ek - m) 



(388) 



Substituting this result into Eqs.(380) and (382), we obtain for the ferromagnetic state 



(70 ( 1 - tanh ( -/3(efe' - /x) 



-^/^y dK'V{k-k') 
Prom this identity equation, we obtain two self-consistent equations as follows 



+ {2S + l)V{Q)n- / dK'V{k-k') 



)nk'. 



lk = \fi j dK'Vik-k'hk'sech" Q/3(efc. - /x) 



(389) 



(390) 
(391) 



Solving Eq.(391) self-consistently, one can obtain the magnitude of the splitting between the up- and down- spin energy 
spectrum. The zeroth order approximation of Eq.(391) yields the celebrated Stoner condition for ferromagnetism, 



1 



=no,/ 

= VF{0) '. 



dK 



dn{ek) 
dek 



VN{0). 



(392) 



Here, F{0), N{0), and V{0) represent the static Lindhard function, the density of states at the Fermi surface and 
the average exchange interaction respectively. The total energy and the thermodynamic potential are obtained in the 
same fashion as for the paramagnetic case. 



E = VJ2 J dK 



nka, 



= --ln2- -y V / dKln 

f3 n ^ I 



l + e 



-/3(efe<T-A') 



where 



^ka = lv{0)n - i / dK'V{k - k')nk'., 



nka 



1 -I- e0(ikcr-ij) ' 



(393) 
(394) 

(395) 
(396) 



Cka and Tikij can be determined self-consistently once 7fe is obtained. 

In this chapter, we have shown that the functional Schrodinger picture formalism can produce the well-known 
Hartree-Fock or mean field results through variational process with Gaussian trial functional. 

We also have shown that, in certain special cases, a shifted Gaussian can be employed to obtain higher terms 
which do not appear in the Gaussian approximation. Although the present method is quite versatile and can be 
applied to any models, including the electron gas, the Hubbard, and the nonlinear sigma model [4], it is quite clear 
that the variational method lacks means for systematic improvement. Therefore, it is not possible to go beyond the 
Hartree-Fock result within the theoretical schemes presented in this chapter. 

In order to overcome this problem, a time-dependent functional Schrodinger picture theory is constructed [5]. It is 
known that infinitesimal time-dependent Hartree-Fock fluctuations around the Hartree-Fock ground state yield the 
random-phase-approximation (RPA) results (p. 292, Gen. Ref. 2). Indeed, it was shown that, in the small oscillation 
regime around the Hartree-Fock solutions, the RPA results can be obtained. The theory was applied to the electron gas 
to yield the electron-hole excitation spectrum, the collective plasma oscillation, and the spin wave energy. However, 
again it should be noted that this line of approach is limited to the RPA results. Often, for systems with strong 
correlations, the RPA results are not satisfactory. Therefore, it is necessary to provide a theoretical scheme which 
can go beyond the Gaussian and the RPA results. 

In the next chapter, we discuss such efforts using the functional Schrodinger picture. 
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IV. FURTHER DEVELOPMENT OF THE FUNCTIONAL SCHRODINGER PICTURE 

METHOD 



A basic problem inherent to any variational calculation is that it lacks means for systematic improvement. 
Therefore, in order to devise a many-body theoretical scheme applicable to wide range of problems, it is necessary to 
overcome this problem. 

In undergraduate quantum mechanics courses, students are taught generally two types of approximation schemes, 
perturbation and variational methods. Each method is known to have its own advantages and limitations. Pertur- 
bation theory is applicable only when the perturbing term is small, whereas the variational method lacks systematic 
means for improvement as mentioned above. Efforts to combine these two approximations and provide a new scheme 
which retains the advantages of both methods have been existed in various field of physics. In this chapter, we re- 
view these efforts and introduce the variational perturbation theory and the optimized perturbation theory in plain 
quantum mechanics languages. Then, we construct the variational perturbation scheme based on the functional 
Schrodinger picture. The theory is applied to the A^^ field theory in order to compare the scheme with other existing 
methods. In the remaining section, the theory will be extended to the optimized perturbation with the A^^ as a 
model system. 



A. Variational and Optimized Perturbation Theories in Quantum Mechanics 
Perturbation Theory 

In order to explain the salient features of perturbation theory, we reproduce the discussion on perturbation theory 
from Negele and Orland (Chapt. 2, Gen. Ref. 2). In the process, we will show that how the new concepts of 
variational perturbation and optimized perturbation are related to plain perturbation theory. 

Perturbation theory is based on the belief that the behavior of a physical system is continuous in some 'small' 
parameter describing the difference between a solvable problem and the actual system. It is crucial to note, however, 
that in general perturbation theory yields an asymptotic rather than a convergent series. Only under appropriate 
circumstances, it may yield a useful physical approximation, but it cannot be applied systematically to arbitrary 
precision. 

The salient feature of asymptotic expansions may be illustrated by the following simple integral 

Z[g) = [ ^e--V2-fl.V4 (397) 
J v27r 

corresponding to the classical partition function Z = J rfxe"^^^^ for a particle in a potential, V{x) = /2 -\- gx^/^. 
Physically, the classical or quantum behavior in this potential changes completely when g changes sign, since a particle 
is localized when <? > and not when <? < 0. Thus, we expect that the theory is nonanalytic at (/ = and that an 
expansion in powers of g has zero radius of convergence. 

The perturbation series for Z{g) is easily obtained by expanding e~^^ with the result 

Z{g) = Y,9''Zn (398) 

n 

with 

^ (-g)"(4n-l)!! 
n!4" 

_ (-g)"(4n)! 
n!l6"(2n)! 

. 1 fl^y, „^oo (399) 

where the asymptotic behavior in the last line is obtained using Stirling's formula n! = •\/27rn"+^e^". Since Z^ grows 
like n\, the series diverges as expected from the non-analyticity at g = 0. 
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The crucial point for the present discussion is the fact that, under appropriate circumstances, a finite number of 
terms of an asymptotic series may give an excellent approximation. The residual error after n terms is bounded by 
the (n + l)st term 



Rn = 



z{g) - Yl 5™^" 



m=0 



/ 
/ 



dx 



'2lT 

dx 



^-x^/2\^-gx^/4 



m.\ [ 4 



m=0 

4 X n+1 



/27r (n + 1)! V 4 

— I 7 I 



(400) 



so the approximation continues improve as long as g"'Zn decreases. Fig. IV. 1 shows i?„ and as a function of 

n for several values of the coupling constant g. From the asymptotic expression for g^Z^, Eq.(399), we note that the 
minimum occurs for n ~ so that the minimum error is ^"^nlmm ^ -\/45'/7re~^/^^. This exponential dependence 
on the inverse coupling constant is characteristics of perturbation theory for weak coupling. 

The real problem comes at large coupling constant where the series begins to diverge in very low order. At 5 = 0.1, 
for example, one must stop at the third term after only attaining an accuracy of R3/R0 = 25%. 



10' 
10° 
10-' 
10-2 
10-3 
10-" 

10-5- 
0-6 

10-M 

10-8 

10-3- 

10-10- 

10-11- 
10-'2_ 

10-13- 



g=o.2 ■ 




g=0.01 



10 



20 

n 



30 



40 



Fig. IV.l Asymptotic expansion of Z{g). For five values of the coupling constant, ^g"\Zn\ is plotted at each integer 

n and the residual error R„ is plotted at each half integer n + i. The arrow symbolizes the variation process carried 
out for the variational perturbation theory, whereas the solid line connecting the minima corresponds to the optimized 
perturbation approximation, (from Gen. Ref. 2). 
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In order to overcome this problem, many prescriptions including the Borel summation (Section 7.5, Gen. Rcf. 2) 
have been proposed. In this lecture, we discuss two systematic approaches to overcome this problem, namely the 
variational perturbation theory and the optimized perturbation theory. Fig. IV. 1, which is Fig. II. 2 in General Refer- 
ence 2, is reproduced here to explain the relationship between the plain perturbation and the two methods mentioned 
above. 

Variational Perturbation Theory 

Here, we explain the concept of variational perturbation theory using a simple quantum mechanical language [1]. 
Although this concept has been studied intensively in connection with path integrals [2], field theories and other 
physics subfields [3] , we believe that it had never been discussed in a plain quantum mechanical language before Ref . 
1. 

In a variational calculation, one first chooses a trial function ^'(A) as a function of a variational parameter A for a 
given Hamiltonian H. Then, the ground state energy is estimated by minimizing the expectation value, (^'aI-H^I^'a) 
against A. 

In the present formalism, we also first choose a trial function \l/n(A), but we require H\\^n{^)) = -Eii(A)|\l/„(A)), 
where H\ is called a parent Hamiltonian. We will discuss implications of this requirement below. The original 
Hamiltonian is now rewritten 



H 



Hx + H- 
Hx + H', 



(401) 



where H' = H — Hx is the new renormalized perturbing Hamiltonian. Clearly, success of the following perturbation 
calculation depends how small H' can be made. In order to obtain an optimum Hx, we first determine A through the 
condition, 



5(M/„(A)|ir|M/„(A)) =0. 



(402) 



Here, we note that this condition is not limited to the ground state only, but is valid for all excited states. Therefore, 
A will be generally dependent on the state number n, of which energy we want to estimate. 
The standard Rayleigh-Schrodinger perturbation expansion is given by [4] 



1=0 L 

E^ = E(^^+^('\n\Hj 



1=0 



Q n 



l—(E^-E^^)-Hi) 



(403) 
(404) 



where Q„ = 1 - P„ 1 - |n)(o)("Vl- 

In the present theoretical scheme, we use Hq 
order is given by 



Hx and H' = H — Hx ■ The energy expression up to the second 



En= {^n{Xn)\HxJ^n{\n)) 
+ {^n{Xn)\H'\^n{Xn)) 

^^|(*.(A„)|F'|*„(A„))|2 

+ ■■■ 



4°^(A„) 



(405) 



where lfA„ |^'n(A„)) = En\x„)\^/ n{Ki)) ■ We note that the first two terms correspond to the variational calcula- 
tion result (4'„(A)|iJ|4'„(A))A=A„ with A„ determined through minimization. The third term is the second order 
perturbation term with a renormalized perturbation which provides a systematic improvement over the simple vari- 
ational ground state energy. The perturbation expansion, Eq.(405) differs from that obtained through conventional 
perturbation theory. The basis function used in the perturbation are those obtained through a variational process. 

As a simple example of the variational perturbation scheme, we consider the anharmonic oscillator problem. We 
will sec that seemingly simple example exhibits all the salient features of the variational perturbation theory. 

The Hamiltonian is given by 
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H=^ + ^x' + bx' , (406) 

2m 2 

where b is positive. An obvious choice for the trial function is given by 

1 

*n(a;) = ( j e , (407) 

for the ground state, where is the variational parameter. The corresponding parent Hamiltonian is given by 

Using this parent Hamiltonian, the original Hamiltonian is rewritten as 

We denote the nth eigenstate of ilo by |nn). Then, the expectation value of (na|if|na) is given by 

(nn|-ff|nn) 
= (nQ|i?n|nn) + (n|if'|n) 

= (nnl-ffjilno) - "^(^ ^ ) + fo(nn|a;''|nn) 

= *<" + 5) - ^^^^<'» + ^' + 4^<'"' ■ 

Here, we used the standard quantum mechanical results [5] 

(nalx^lna) = (^^^ (Gn^ + Gn + S) . (411) 
Taking a variation on (nn|il|nn), we obtain a relation which determines f2„, 

O3_.2^,„-___^=0 . (412) 



TABLE I. Comparison of the ground state energies(in eV) obtained using different approximation schemes. The values in 
the parentheses are the ratios to the exact values. The values of imf2§ are also shown as references. 



b(eVA-^) 





0.01 


0.05 


0.25 


Perturbation 


1.4318427 


1.5279252 


does not converge. 


Theory 


(99.935%) 


(95.962%) 




Variational 


1.4333279 


1.5968858 


2.0664772 


Calculation 


(100.038%) 


(100.293%) 


(100.929%) 


Present 


1.4327276 


1.5912088 


2.0412648 


Method 


(99.997%) 


(99.937%) 


(99.679%) 


Exact 


1.4327725 


1.5922195 


2.0474629 


Energy Value 










0.5770839 


0.8227827 


1.6423320 
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Substituting these results into Eq.(405), we obtain 

AvrfiVLl 



En = ^(2n+ 1) - —-i2n' + 2n + l) 



1 I bfi^ \ 64n5 + leOn^ - 336n3 - 664n2 - 280n - 24 

\Aw?ni) (2n + l)2 ■ ^ ' 



<Problem> Show that En to the second order is given by Eq.(413). 



In order to show the systematic improvement achieved by the variational perturbation theory over the variational 

2 

or the plain perturbation theory, we have carried out numerical calculations. For this purpose, we choose ^^^y" = 
0.beVA~'^ and carry out calculations for various values of h. The results for the ground state energy are given in Table 
4.1 [Table 1, Ref.l]. The result shows that the present method is clearly superior to the conventional perturbation 
theory and provides a systematic mean for improvement over the variational calculation. First of all, the present 
method gives highly accurate values in the regime where the conventional perturbation theory is not applicable. This 
is because the pertm-bing Hamiltonian has been renormalized through the variational process. 

Before we discuss the convergence problem of the variational perturbation theory, we briefly mention that the 
method can be applied problems other than anharmonic potentials such as helium atom and also to excited states 
[1]. Also, any well-behaving potential can be expanded around local minima 

V{x) = V{xo) + ;^^^"H^)U=xo(^ - ^o)". (414) 

n=2 

Therefore, we observe that above discussion on anharmonic oscillator is generally applicable to any well-behaving 
potentials. In the next section, we will show that this approximation corresponds to the variational perturbation 
process based on the Gaussian basis. 

Next, we discuss the convergence problem of the variational perturbation process. Recently, this problem has been 
studied in detail and has been shown that the expansion does not converge at higher orders for strong coupling 
constants [6]. In fact, this behavior is not totally unexpected. 

In the early part of this section, we have shown that perturbation theory is only asymtotically convergent and even- 
tually divergent. The variational perturbation is also basically a perturbation expansion only with a renormalized 
perturbation. This process is explained on Fig.IV.l. As indicated by the arrow, the variational process renormalizes 
the the coupling constant first. Then, perturbation calculation is carried out with the renormalized coupling constant. 
As shown by the anharmonic example, this renormalized perturbation generally allows one to obtain much higher 
accuracy than allowed in the original perturbation expansion, although the variational perturbation theory will also 
eventually diverges at higher orders. 

Optimized Perturbation Theory 

Optimized perturbation theory was proposed to overcome the convergence problem in the perturbation theory [7, 8]. 
It may look formidable in the field theory, but, in quantum mechanics, it is a rather simple and straightforward idea. 
In the optimized perturbation expansion, the variational process is not made in the beginning, but at the last stage 
of the perturbation calculation. In the process, the Hamiltonian is divided as in the variational perturbation method. 

And the perturbation calculation is carried out using Eqs.(403) and (404). Variation is made on E'-'HX), the Ith. order 

perturbation result, and, then the obtained A„ is substituted into En\x). Therefore, the value of An' is different on 
different I. 

We believe the method is best understood through a simple example of anharmonic oscillator [9]. 
We use the same Hamiltonian, Eq.(409) 

H = Hn+H' 
and the same trial function, Eq.(407) 
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We carry out the perturbation calculation on En using Eq.(404) but without the variational condition, Eq(412). Then, 
to the second order, we obtain 

En = {na\Ha\nn) + {na\H'\na) + ^ '^^J^^oT 

fcytn En Ef. 

= ^(2n+l)-^^^(2n + l) 

+4;;^(^" 1603 (^--+1) 

-g^^^(34n3 + 51n2 + 59n + 21). (415) 

We find that, by taking a variation on En with respect to O, dEn/dVt = 0, O does not any real roots. In order to 
cope with this problem, we examine the Hamiltonian, Eq.(409) and the perturbation expansion, Eq.(404). We note 
that, if the perturbation calculation is to be carried out to infinite order, then, the final result will be independent of 
n, since the original Hamiltonian does not have 0. The O dependence in the energy values comes from 'incomplete' 
perturbation calculation. 

Therefore, unlike in the variational perturbation where one looks for minimum, here, we search for a point where the 
energy value is least sensitive to the variational parameter O. This is called 'the principle of minimal sensitivity' [8,9]. 
In the principle of minimal sensitivity, one takes the point which satisfies d^En/dO? = 0, as the point of the least 
sensitivity, when there is no real value to satisfy dEn/ d^l = 0. It will be shown that in the odd order approximation 
, there exists at least one real value of fi, whereas in the even order, one should choose an O value which satisfies 
d'^En/dfl'^ = 0. This comes from the fact that, at the odd order, the highest order of O is even and odd for the even 
order. 

Therefore we take d^En/dfl^ = for Eq.(415) and, by numerically solving the equation, obtain flig. The value of 
Clis is substituted into Eq.(415) to obtain the energy value £„. Here, wc note that, in the optimized perturbation 
theory, it is not possible to obtain an analytic expression for the energy unlike in the variational perturbation theory, 
Eq.(413). 

In order to study the convergence problem of the three perturbation theories, we have carried out the optimized 
perturbation calculation up to the third order, 

En = {nn\Hn\nn) + {n\H'\n) 
\{kn\H'\nn)\' 



E 



{na\H'\kn){ka\H'\l^){ln\H'\n\n) 

^1^1^ (jpiO) _ p,(0) ^ , p(0) _ p(0) X 

E \{ht\ir\iiii)\-{iu_i\ir\iiii) 

kjin - J^k ) 

= -(2n+l)-^^(2n+l) 



(34n^ + 51n^ + 59n + 21) 



- 320^ + 1) + 4rn2^a + 6^ + 3) 
16rn407 (170n^ + 255n" + 292n + 105) 
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+ ^(375n* + 750n^ + 1416n^ + 1041n + 333). (416) 

l6m°U° 

<Problem> Prove Eq.(415) and (416) . 

This equation has a real solution for dEn/dQ = 0, which can be substituted into the equation to obtain the energy 
value. 

In Fig. IV. 2, we compare the three perturbation methods for the ground state(n = 0) with parameters = 
0.5eVA^^ and 6 = 0.01eyA~'*.Thc figure wliows clearly that the optimized perturbation has the best convergence. 
Now, we consider the relation between the optimized perturbation theory and the conventional perturbation. 





t 1.437364 (100.32%) 


• perturbation theory 

o V-P method 

T Optimized expansion 




•.1.433278 (100.035%) 


• 1.433051 (100.019%) 


Exact 




01-432782 (100.0007%) 
. ; . T 1 .432773 (1 00.00005%) 


value 






: 1.432772 
(100%) 


;t1.. 432769 (99.9998%) 
■ . 1 .432727 (99.997%) 

• 1.431842 (99.935%) 




1 1 
1 2 


1 1 
3 4 



Fig. IV. 2 The three different perturbation theories are compared for the parameters mui /2 = O.SeVA-^ and b = 
O.QleVA~* up to the third order. 



The optimized perturbation theory renormalizes the perturbation(coupling constant) at each order and choose the 
minimum(or near minimum) point for E^. Therefore, the optimized perturbation theory can be pictorially represented 
by solid line in Fig.IV.l. 

Here, it should be mentioned that although the optimized perturbation theory may solve the problem of convergence, 
its practical use is rather limited due to the complicated nature of the detailed calculation. Quite often, the one-time 
renormalization scheme of the variational perturbation theory allows enough accuracy required. 

The variational and the optimized perturbation theories in quantum mechanics via path integral have been inten- 
sively studied by Kleinert and his colleagues and any higher terms can be readily generated through computer [2]. 

In the remaining sections, the variation and the optimized perturbation theories will be applied to bosonic field 
theories using the A(^^ model as a model system. 



B. Variational Perturbation Theory in The Functional Schrodinger Picture 

It has been shown that in Chapter III, Gaussian approximations invariably lead into mean-field results. In order 
to go beyond Gaussian results, various approximations schemes including functional integrals and background field 
methods have been proposed[l'~3]. Among these methods, except several variational calculations with correlated or 
non-Gaussian trial functions, majority of the schemes are basically based on the concepts of the variational or the 
optimized perturbation theory. Therefore, here again the approaches going beyond the Gaussian approximations can 
be classified by the methodologies employed, namely, the operator (Green's function method), the functional integrals, 
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and the functional Schrodingcr picture methods. It has been observed that the first two methods naturally lead into 
Feynman diagrams again albeit renormalized and somewhat simplified. Therefore, in order to go beyond the second 
order, careful enumerations on inequivalent diagrams are still required. However, functional Schrodinger picture does 
not require such counting processes. Being a basic perturbtion theory but with a Gaussian basis, it is enough to carry 
out the standard Rayleigh-Schrodinger perturbation calculation given by Eqs.(403) and (404). The only modification 
is that the basis are now functional instead of functions. 

Here, we present a bosonic Rayleigh-Schrodinger perturbation theory based on Gaussian wavefunctional. We 
consider a model with the Lagrangian density 

^ = Id^^.d'^^x - ^('/>.), (417) 

where x = (x^,a;^,- • •,.t^) represents a position in D-dimensional space and (j>x = (t>{x) is the field at x. We note 
that this model can be readily reduced to non-relativistic many-body theory by changing the kinetic part. Also, the 
theory can be reduced to non-relativistic quantum mechanics by requiring _D = 0. 
For the variational perturbation process, we use the Gaussian trial functional 



|0) 



(0) 



= TVexp |i j Vy(j)y J {4>x- v)fxy{(t>y - V)], (418) 

where N is the normalized constant,/^ y ^ -f '^^'"^ '^^V^ ^^'^ f^v ~ J d^pf{^)e^^^'^~^^ with p = {p^jp"^,- ■ -,p^)- The 
classical constant is equal to the Gaussian vacuum expectation value of 

<p= W(0|</.,|0)(0). (419) 

The Hamiltonian density corresponding to Eq.(417) is given by 

W = TTxdtCpx - 

= l7Tl + ^{dx<t>xf + V{ct>x), (420) 

where tt^; = as shown in Eq.(190). Here, we set c = ?i = 1 following the field theoretical notation. In order to 
calculate the ground energy state density e = {H), it is necessary to evaluate the following matrix elements[4], 

(°)(0|(a,<^,)2|0)W = {dx<pf - {{<f>x - V)dl{4>x - f)) 

^2 , 1 /• d^P P" 



where we utilized the Gaussian functional integrals, Eq.(154). We also have 

■°'»"i»>'°' = ((i^)'> 



I 

,2,1/' d'^p p2 



('')(0|F(*,)|0)(») = r * 

J-oo V^TT 

J-oo v27r 



oo 



^V(3) r _^^e-«'/2e^/^(-T^^^ '^e'^'f 
J-oo v27r J-oo v27r 

da -^v2t./ a rzi \ 



/oo T I 
^e-""v(a^fxx'+<p). (423) 
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<Problem> Prove Eqs.(421)-(423). 



Collecting the results, we obtain for the ground state energy density, which is also called effective potential in the 
field theory, for the Gaussian trial functional. 

e|v,P„/| s <"'(0[«,|0)<»' 



2V ■ 2 ^ ■ 4 

1 f d'^p p2 



+ 



4 J (27r)^/(p) 



■if 



(424) 



As the first process of the variational perturbation theory, we take variations on e with respect to Vx and f{p). 
We observe that 6e/5p = simply gives = as expected. 6e/6f{p) = gives the crucial minimizing condition 
corresponding to Eq.(412) in quantum mechanics. 



5e 



= = 



+ 



+ 



P 



1 



4(27r)^ 4(27r)^ V P{p) 
1 



df5V{(i)e 



(27r^) \p{p) 



(425) 



where we used the second relation in Eq.(423). This equation can be arranged as 



\p)=p'' + J 

J-oo 



2 , 2 
= P +/" , 



(426) 



where y(")(0) = d''V{z)/dz" = J^^{dq/V2TT){iq)"V{q)e'i\ 
The renormalized mass or coupling constant fi^ is given by 

Here, it is convenient to introduce the following notation 



> J (27r)^ 



(p2 + Q2)n 



(427) 



(428) 



Note that f'^ = hin"^). 

Substituting the above relations, Eq(426)'~(428), into Eq.(424) and noting that d^/dx = 0, the Gaussian effective 
potential becomes 



Vg(^) = ^Io(m')-^/i(m') 



+ 



(429) 



We note that when Vg{v) has the absolute minimum at (po, f^ifo) becomes the physical mass(or the renormalized 
mass) and Vcif) represents the vacuum(ground) state energy. The symmetry of the ground state can be discussed 
using Eq.(429). 

For the second stage of the variational perturbation theory, it is necessary to partition the Hamiltonian in two 
parts as we did in Eq.(401) and (409). This can be best achieved constructing annihilation and creation operators for 
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Gaussian particles (Gaussions), of which the Gaussian functions are exact solutions. 
Following Ref.[4], we define 



Af{p) = 
A\{p) = 



1 



2(27r)^/(p) 
1 

2(27r)^/(p) 



1/2 



1/2 



e-^P-[f{p){<P,-^)+in.. 
I e'^'^[/(p)(</..-<^)-m,] 

J X 



(430) 



We can readily show that these operators satisfy [Af{p),A^^{p')] = 6{p' — p). Based on these operators, one can 
construct the parent Hamiltonian 



Ha 



dp f{p)A\{p)Af{p) 



(431) 



<Problem> Show that [Af{p),A\{p')] = 5{p' -p) and prove Eq.(431). 



The Gaussian ground state(vacuuin), Eq.(418) is the ground state of Hq with the zero energy eigenvalue e'q \ The 
exited states of Hq are 



1 " 

\n)^''^ = ^MA]{Pi)W\ n=l,2,. 
with the corresponding energy eigenvalues 

n 
i=l 

jn)^"^ are normalized, 

■ Pi(n) k=l 



(432) 



(433) 



(434) 



where the Pi{n) represents a permutation of the set {«/j}={l,2, ■ ■ -,71] and the summation is over all Pj(n)s. \nY^^ 
describes a n-particle state with the continuous momenta pi,p2,---,Pn- \0)'^^^ and | n) with n = 1 , 2 , • • • , 00 constitute 
the complete set for Hq- 

We note that the above Gaussian formalism is quite general and convenient for any system with a well behaving 
potential V{(l)x)- Now, the Hamiltonian is rewritten H = Hq + Hj = Hq + {H — Hq) with 



Hi 



Following the Rayleigh-Schrodinger perturbation procedures, Eq.(403) and (404), we write 



1=0 L 



(£;„ - - Hi) 



In) 



(0) 



and 



1=0 



Ho - El^> 



In) 



(0) 



(435) 



(436) 



(437) 



where Q„ = / d^Pid^P2 • • • d^pj\j)(^^(^Hj\. 

For the case, n = and I = 0, Eq.(437) gives the vacuum(ground state) energy up the first order Eq = Eq^^ + Eq'^ 
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fx VG(iy5), which is just the product of the Gaussian effective potential and the space volume. 

Application to the A^^ Field Theory. 

As an example, we consider the potential, V{4>x) ~ \^'^4''i + (''^/4!)0x) which was widely studied in connection with 
the Gaussian approximation [5]. Thus, this model allows direct comparisons with existing results and can be readily 
reduced to non-relativistic many-body theory and quantum mechanics. 
With the above potential, Eq.(427) becomes 



= m2 + iA^2^iA/i(M'), (438) 
where we used ^a^^+^e""' = and ^a2"e-«' = 2"" • 1 • 3 • 5 • • • (2n - 1). Also, Eq.(429) becomes 

Vg{^) = ^mV' + ^A^4 + ijo(M) - 3^A/f (m^), (439) 

which is the standard Gaussian effective potential for the Xcj)^ model. 

<Problem> Prove Eqs.(438) and (439). 

To obtain the effective potential of the A^^ model up to a given order using Eqs.(436) and (437), we need the 
following matrix elements: 

(°)(n|i?.|n)(°)=nO|i?,|0)(°)(°V|n)(°^ + ;^T7^ E ^^(p^-^) 

+P'j^-PHn-l) -Pin) 

><[fiPL.JfiP'jJfiPH.-r,)fiPiJ]-'^ . (440) 



Pij(n-2) fc=l 
(n-i)^ 



X E riA(K,-4)c5(p,J[/(p,J]-ij, (441) 

^^^'kn-^) -Pin-2 -PHu-I) ' Pin) 

X[/(Pi(n-.))/(P^n-)/(P»(n-l,)/(PiJ]"^ , (442) 

v/^^K^^ 2^2(2^, ^^|;^_3^i^^^ 

x[/(Pi(„_„)/(Pi(„_„)/(ftJ]-^ (443) 
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>«^Kn-3) +Pin-2 + + Pin) 

X[/(Pi(„_3,)/(Pi(„_.))/(Pi(„_,))/feJ]-^ , (444) 

with 



for fc < 

MP^.-P'J={ 1, for k = 

for fc > 




Here, the index ik G {1, 2, • • • , n} with fc = 1, 2, • • • , n corresponds to |n/°^ and jk S {1, 2, • • • , n'} with = 1, 2, • • • , n' 

to *^''-l(n'|. Pij{l) represents a given permutation of I momenta pi-^,pi^, ■ ■ ■ ,pi^ paired respectively with Pj-^tPj.^, ' ' ' jPjji 
and X]p^^.(;) is over all different Pij(l)s. For any Pij{l), H,«2,''')*( are different from one another, and so are 

<Problem> Prove Eqs.(440)~(444). For the calculation, it is much more convenient to carry out the normal 
ordering process beforehand [See Appendix C.]. 

Employing the above matrix elements, a straightforward, yet lengthy calculation according to Eq.(437) gives the 
effective potential of the A^^ field theory up to the third order as 



with 



(277)^5(0) 

•>i/\ 1-'- 2/2 ^ 1\'2 ^ ^ 2 ^ ^ 
+ -3^) 

+ 128 ,r-sn ^+^2^ ^'''^ 



^ = / (2fF (4iT[/i(x)/i(y)/i(x + W + fi{y) + h{x + , 

^1 = / (2^ (2^ [A + y)]-'[^ + h {x) + h {y) + h{x + y)]-\ 

^2 = / (^iT (4TT[/i(x)/i(?y)/i(x + y)]-'[h{x) + My) + h{x + y)\-' 

x[l + /i(x) + /i(y) + /i(,x + y)]-i , 
B = S i0^j0^j0^[h{x)h{y)h{z)h{x + y + 

X [h{x) + + h{z) + h{x + y + z)]-^ , 

^1 = / i0^^j0^[hi.x)h{y)h{z)h{x + + y + ^)]-' 

X [h{x) + + + y)]-^[h{x) + + h{z) + h{x + y + z)]-^ 
B2 = I i0^^j0^[h{x)h{y)h{z)h{x + y)h{x + 

X [/iW + + + + /iW + + 

C = !j0^j0^^^[h{x)h{y)h{z)h{u:)h{x + y + z)h{x + y + u:)]-^ 

X [/i(x) + + h{z) + .h{x + y + z)]-^[,h{x) + h{y) + /i(c^) + /i(.t + ?y + co)]-i 



and = vl + u^. In Eq.(445), the second, third, and fourth terms are the second order corrections and the 

last three terms represent the third order corrections to the Gaussian effective potential. 

<Problem.> Prove Eq.(445). 

Quite surprisingly, we note that the second order result of Eq.(445) has an additional term, — [/i^ — 

(A/3)(p^]^ which does not appear in other references [6]. In order to see the origin of this term, we note that the 
normal ordered form of Hi is given by [See Eq.(513).] 

+ ^^:{<i>.-^f--+^,:{r-vf:], (446) 
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where : : denotes normal ordering. The additional term arises from the second term, which disappears if we choose 
< (f'x >= However, we note that this result is from the simple Gaussian calculation and does not include higher 
order contribution, therefore, when higher order contribution is included self-consistently, it is expected that this 
term produces nonvanishing contribution. 

It should be noted that in obtaining the above results, no diagrammic considerations are made. Although the 
calculation maybe lengthy, it is rather straightforward and conceptually simple. We believe that the complexity in 
selecting relevant diagrams and corresponding symmetry factors have, so far, prevented obtaining the third order 
results using the operator and the function integral method based on the Heisenborg picture. 

Before concluding this section, we compare the present results with the ones in the quantum mechanics, Eqs.(413) 
and (416). Quantum mechanics requires D = ; ip{x, y) x{t), d^p = 1 and {2n)^ = 1. Also, we observe = for 
anharmonic oscillator. Then, Eq.(445) is reduced to become 

VG(0) = i/o(/x2)-3^A/2(/.2) 

384 ^ 512 ii» 
_ 1 _ 3 6 _ 362 27 6^ 

Here, we used /q = = = 3'^ ~ TS ^^'^ ^ ~ A/46. Eq.(447) can be directly compared with Eq.(413) and 
(416) and confirms the correctness of the present calculation. 

The discussion in this section has been confined to the bosonic case. In order that the variational perturbation 
theory can be applied to condensed matter systems, it is necessary to develop Rayleigh-Schrodinger perturbation 
theory for fermionic functionals, which is a future subject. 



C. Optimized Pertubation Theory in The Functional Schrodinger Picture 

In this section, wc present an optimized perturbation theory in the functional Schrodinger picture with an external 
source[l]. First, a free field theory with an external source will be constructed. Based on this free-field basis, a general 
expression on the optimized effective potential for a class of model potentials will be given[l,2]. Then, the theory will 
be apphed to the X(j)^ model. 

The Hamiltonian with an external source, Jx = J{x), is given by 



■(t>x - \fxx + \ I Jxh-^Jy\ , (448) 



where f^y = {\/-dl + fi'^)S{x-y) with fxzfzy = ^ix-y) and h^y = {-dl + ii^)5{x-y) with Kzh^y = 5{x-y). 
The last term is introduced to make £'o[J] = in the presence of the external source term —Ja;(/)a; (Compare with 
Eq.(429).). The functional Schrodinger equation is readily solved using the power-counting method. The ground state 
wavefunctional is a Gaussian-type functional 

|0; J)W = TVexp { - i ^ ^ (^c/., - ^ /i" V,) - ^ h-^J^ }, (449) 

where is the normahzation constant. It should be noted that, here, (°)(J; 0|</)a;|0; J)^"^ = h~lJz which, unlike in 
Section IV. B, is dependent on x. Then, accordingly, annihilation and creation operators can be constructed as 

= ( 2(27r)W(p) ) i ""'''[■^(^')('^- - / "^^) + ^'^-l' (^50) 

and 
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with [Af{p; J), A\{p/; J)] = 5{pi - p) and Af{p- J)|0)(°) = 0. 

<Problem.> Show that H^' ^ = J d'^p f{p)A\{p; J)Af{p; J), where f{p) = ^p^ + m^- 
The eigenwavefunctionals for exited states can be written as 

1 " 

|n; J) W = ^ n •^)^°^' " = 1' 2' • 



, oo 



(452) 



and the corresponding eingenergies are 



(453) 



The eigenwavefunctionals |n;J)'^°) and |0;J)'-°-' arc orthogonal and normalized, ^'^^ ( J; m|n; J)*^"-* = 
<5mn;;7[ (n) 11^=1 "^(Pfe ~Pik)- Hcrc Pi{n) represents a permutation of the set {ik] = {1, 2, • • •, n} and the summation 
is over all Pi(n)'s. \n; J)^"^ describe a n-particle state with the continuous momenta pi,p2, ■■■,Pn- |0; J)^°' and \n; J)^°^ 
with n = 1, 2, • • •, oo constitute the complete set for Hq ^ and satisfy the closure |0; J)^°^(°^(J; 0| +^^;^ J d^pid^p2 ■ 
• •dV>;J)(")(")(J;n| = 1. 

We now introduce the model potential through the following scalar Lagrangian density 



(454) 



We assume that V{(l)x) = J -^=V{Q)e^^^'^'-' , at least, in a sense of tempered distributions [3]. In fact, this assumption 
also applies to the previous section. For the system, Eq.(454). the time-independent functional Schrodinger equation 
in the presence of an external source of Jx is 



H 



) = En[J]\^n) 



with the Hamiltonian H = /^[^tt^ + \ {dx4>l) + 

Here, the eigenvalue -En [J] is a functional of Jx- For convenience in the calculation, we shift (jix 
constant). Thus, the above equation is modified 



H{cj>x + $) 



J X 



+ $) 



*„[(/.^ + $, J] 



= .E„[ J; + J]. 



(455) 
$($ is a 

(456) 



This shift is made in the spirit of the background field method[4]. However, we will not go into the background field 
method. 

Therefore, the above shift can be regarded as a just mathematical convenience. Normal ordering the Hamiltonian 

with respect to a normal-ordering mass M (Sec Appendix C.) and inserting a vanishing term /^[l/U' 
jjL an arbitrary mass parameter into the Hamiltonian, we have 



i-^-^^-iMV^lwith 



M 



with 



and 



J,fX 



H{cj>x + ^)- I J.((/'. + $) 

J X 

= H^'^ + iff* - C 



.i^,^cl>l + NM[V{ct>x + m 



} 



C = J^ -lUx + lJ^ Jxh-^Jy + lloiM') - ^h{M') + Jx^ 



(457) 



(458) 



(459) 
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Here, the normal ordering of the shifted potential is given by[5] 



27r 



(460) 



For IniM"^), sec Eq.(426). Since C is a constant, we can now carry out the Raylcigh-Schrodingcr perturbation 
calculation only for iJj To check the order of the perturbation calculation, an index factor 5 will be attached to 
Hj' , so that Eq.(456) is modified to become 



[ffo''^ + 5fff*]*„[0, + a>,J;<5] 
(£:„[J;$,5]+C)*„[<^, + $,J; S]. 



(461) 



Now, we carry out the optimized perturbation expansion as explained in Section IV. A. Obviously, the zeroth-order 
approximation to Eo[J; $, 5], i?o"^ [J; ^] satisfies 



eI°\j,^ + c = e^- 



(0) 



0. 



(462) 



To the nth order of 6, 

4°) [J; $]=(«)( J; 0|JJf * 



Oo- 



(0)r 



|0;J)(°), 



(463) 



with Q„ = „ / rf^Pid^P2 • • • d'^pjlj; J)W W(J; j|. 

Thus E,[J;<S>,5]=E^"\j;<P]+Y.n=iS''4"\j-^'^]- 

Before carrying out the perturbation calculation, we note that a Legendre transformation of Eq[J; $, 6] yields the 

static eff'ective action[6] 



r4^;$,5] = -i;o[J;$,<5]- / J^v'x. 



(464) 



To calculate the eff'ective potential, one can conveniently take 
Then, the effective potential is given by 



$ to fix the arbitrary shifted parameter $[4]. 



V($) 



L 



(465) 



(^=^,(5=1 



When truncated at a given order of 5, V($) will depend on fj,. To obtain an approximated effective potential, one 

determines /i according to the principle of minimal sensitivity explained in Section IV. A. 

The matrix elements, which appear in Eq.(463), involve only Gaussian integrals except commutators of creation 
and annihilation operators and, thus, can be readily calculated as follows, 

vn!m! rr: 



n m n—i 

j=n— j+l fc=n— i+1 1=1 



(466) 



with n < m. In Eq.(466), V''{z) = ^j^^- For simplicity, in getting the above results, we have employed the 

permutation symmetry of momenta in Eq.(463) for various products of 5 functions. Note that matrix elements of 0^. 
are special cases of Eq.(466). Substituting the above matrix elements into Eq.(463), one can obtain the first- and the 
second-order corrections to Eq [J; (f)] as 
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+ r ^e-"V(a^//.-i-/i(M2) + $ 



(467) 



and 



^ i!2i 7 (27r)^-^ ni=i /(P;fc)(ELi /(Pfc)) 



(468) 



respectively. Here, "| • • • |" represents the absolute value. 

Next, we extract the approximated effective potential for the system order by order. 

At the zeroth order, E^Q\j;(t)] = -C from Eq.(462), and so, taking - ^-^^j^ = JyKyJy + ^ = "P^^ as one 
has J(°' = 0. Consequently, the effective potential at the zeroth-order of 5 is 



V(o)($) 



ri°)[^;$,5] 



fx 



1 1 ]lf2 

-U.--Io{M') + —h{M'). 



(469) 



Up to the first order (Hereafter, any Greek-number superscript, such as "/", "J/", means "up to the order whose 
number is consistent with the Greek number".), 



(470) 



and -^%^=^^ = $ yields 



Jy Jy y Jz 

- ;^e-«V(^) (aV/.-.i-7i(M2) + $ + ^ /^- V.) }= . 



(471) 



When extracting the effective potential up to first order, Eqs.(464), (465) and (470) imply that only the up to 
the first order, J^, is necessary. Owing to j'"^ = 0, it is enough to take Jx = ^ for the last term in the left hand of 
Eq.(471). Thus, can be solved from Eq.(471) as 



(472) 



Even will not be needed to get the effective potential up to the first order of 8, because there exists no linear, but 
the quadratic term of Jx in the zeroth-order term of (£'q[J;$,(5] — J, Jx^) as shown in Eq.(459). In fact, to obtain 



67 



the effective potential up to the nth order, one need the approximated J only up to the (n — l)th order. Now, one 
can write down the effective potential up to the first order 



(473) 



Obviously, this result will yield nothing but the Gaussian effective potential[5]. 

Finally, we consider the second order, ip^ = ip^J - 
however, it is enough to use only for the effective potential. Substituting into Eq.(465), we obtain the effective 
potential for the system, Eq.(454), up to the second order as 



[J'^'^] — $ gg^jj solved for J^^. In the present case. 



+S 



da 



— DC' 

2 



16 



"V(a^/iii-/i(M2) 

' 2/_^^e-«V(^)(aV/.-.^-/,(M^) + $) 

1 



(27r)^ Pip) 



^J\2^ J (27r)0-i)^ 



{f{n-JiPk)+n-Jif{Pk)) 



/(Ei=iPfe)ni=i/fe) 



(474) 



where, one should take 5 = 1 after renormalizing V^^($, 6), and fi is determined from the stationary condition 



= . 



(475) 



Here, V {^) is the effective potential after V ($, 5) is renormalized. If Eq.(475) has no real solutions, /z can be fixed 

by ^^^^ = as explained in Section IV.A. Note that in (1+1) dimensions, {^[P^-IoiM^)] + jM'^IiiM'^)-lii^f-J^} 

and [f~^ — Ii{]VP)] in Eq.(474) with (5 = 1 is finite and, thus, for any (l+l)-dimensional theories which make the 
series in Eq.(474) finite, no renormalization procedure is needed. Similarly, employing Eq.(466), one can obtain higher 
order corrections to the Gaussian effective potential from Eq.(463). 

As an example of the above optimized perturbation theory, we consider the potential, 

F(0,) = ^m^^l + A</.^ , (476) 

which was also studied with the variational perturbation technique. 

Substituting Eq.(476) into Eq.(474), and noting that a2"+ie-«' ^ = for n = 1, 2, • • • and a^^e""" ^ = 

2~" • 1 • 3 • • • (2n — 1) for n = 0, 1, 2, • • •, one can easily obtain the effective potential for the system up to the second 
order 

V''{^,5) = !(/,, - 7o(M2)) + 1m2/i(M2) + <5(lm2$2 + - l^Vx".') 
+5i(/-i - hiM^))[m' + 12A$^ + 3A(/-i - h{M'))] 



/ 



1 



)+fiP2)+f{pi+P2) f{Pl)f{P2)f{Pl +P2) 
d^'pid^p^d^'ps 1 1 



(27r) 



/(Efe=i Pk) + Efe=i f{pk) /(Efe=i Pk) rife^i f{pk) 



(477) 
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Discarding terms with /„(M^)(n = 0,1), the above rcsuh becomes identical to Eq.(2.36) in Ref.7(1990) in Section 
IV. A. This can be verified by carrying out integrations of /„(f2)(n = 0,1) and /''"^(0)(n = 2,3,4) in the above 
reference over one component of each Euclidean momentum. Using jj, fixed at the Gaussian approximation result for 
each order will simultaneously imply that up to each order, the vacuum expectation value of the field operator (f)^ is 
identical to that in the Gaussian approximation. If one chooses to do so, Eq.(477) will yield the second order result 
of the variational perturbation theory in Eq.(443). 

Finally, we point out that taking /x = m will lead to the conventional perturbation result on the effective potential. 



V. SUMMARY AND FUTURE WORKS 



In this lecture note, we have presented a theoretical schemes for study of many-particle systems based on the 
functional Schrodinger picture. In Chapter I, the basic conventional tools for many-body theory including the second 
quantization and the Green's function method are briefly reviewed. In Chapter II, basic ingredients for functional 
Schrodinger picture are treated. They include functional calculus, Grassmann algebra, and the basic functional 
Schodinger picture formalism. In Chapter III, functional Schrodinger picture formalism is applied to various condensed 
matter problems including electron gas, the Hubbard model, the BCS super conductivity, and dilute bose gas with the 
Bose-Einstein condensation. However, the formalism in the Chapter III is limited to the variational approximation 
and shown that it only leads to mean-field(Gaussian) results. 

In Chapter IV, recent efforts to go beyond the Gaussian approximation are presented. Two methods, the variational 
and the optimized perturbation methods are first introduced with quantum mechanical languages. Then, they are 
discussed in depth using scalar field theory. As results of the study, the effective potential (the ground state energy) 
are given, up to the third order in the variational perturbation method and up to the second order in the optimized 
perturbation theory. 

However, it should be noted that no concrete condensed matter models are discussed using above two methods. 

This is because the theory presented in Chapter IV is limited to bosonic case and, thus, cannot be applied to 
electronic systems such as electron gas and the Hubbard model. The work extending the variational and the optimized 
perturbation theories to fermionic system is currently under study. Only when this extension is successful, the 
functional Schrodinger picture approach to many-particle systems can be claimed truly useful. 



VI. APPENDICES 



A. Anticommutation Relations for Fermionic Field Theory 

In order to obtain anticommutation relation for fermion field, it is necessary to start from the spinor field theory 
[1-3]. In general, the fermionic field theory can be studied through either non-hermitian or hermitian field operators. 
Here, we start with a non-hermitian field operators. It is an easy matter to see that the Dirac equation 

[il'^d^ - m)tp = (478) 

follows from the Euler-Lagrange equation 
if we choose 

C = itpj^d^ip — miptp 

= \ [^7''(5^^) - - mi^^, (480) 

where 

tjj = (481) 

and 

A^^'B = ^[A^^'B - {d^A)B]. (482) 
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The 4x4 Dirac matrices are defined by the relation 

= '^9^", (483) 

where the metric tensor gi^" is given by 

/ 1 

9''' =9,^=\l ~0 -1 15 I ' (484) 
\0 -1 

in Minkowski space. 

From C, we find the canonical momentum field 7r(a;) 

nix) = ^=ii^\x). (485) 

This relation naturally leads into the anticommutation relations [1], 

{^P,ix, i), Vl(a^, t)} = S^{x-x')6,j, 

{V'i(f,t), Vi(a?',t)} = = {4{x,t),ijl{x',t)}. (486) 

Note that the above non-hermitian field operators -0 and -0^ are eigenstates of the charge operator. 

Fermionic field theory can also be studied using hermitian field operators [2,3], •4'^{x) = ipix)- Non-hermitian field 
operator can be constructed from the hermitian operators as follows 

'^charge cigenstate V'hcrmitian + ^hermitian- (487) 

In such a case, the Lagrangian density is expressed as 



£ = l^a'^a^V - ^mV/3V, (488) 



where {a\a^} = 26'^ , {a,/3} = O,/?^ = l,a° = 1. 
In this case, we can readily show that 

{Va(x), Vfe(x')}^o=^o' = 5abS{x- f ). (489) 

However, we note that the second half of Eq.(486) is not satisfied in this case. 



References 

[1 ] L. H. Ryder, Quantum Field Theory, Cambridge University Press (London) 1985. 
[2 ] P. Ramond, Field Theory; A Modern Primer, 2nd ed., Addison- Wesley 1989. 
[3 ] J. Schwinger, MIT Lecture Note 1960. 

[4 ] C. Itzykson and J. E. Zuber, Quantum Field Theory, McGraw Hill 1980. 

B. Dual for Gaussian Functionals 

Here, we prove that the dual to a Gaussian functional is given by another Gaussian functional. Specifically we 
prove Eqs.(226)~ (230). 

For a real 6{x) basis, Eq.(223) show that the dual ^[9] of ^[0] is given by 
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ii>[e] = J DO'^^ieY'^, (490) 

where ^'[0] = e^^^^^ . First, the kernel M should have the antisymmetric property to have nonvanishing contribution, 

M* = -M. (491) 

This property can be proved easily as follows, 

eiMij9j = ei{Mj,Y9j = -e,{M^i)'ei 
= -ej{M')jiei = -0i{M%0j. 

Thus, wc have M* = -M . 

Upon substituting -^[9] = ei^^^ into Eq.(490), we have 

J DO'e^^'^'^'e^'^. (492) 



Let 6^ = Ui + KijOj, then 



Thus, 



0'9 = {ui + Kij9j)9i = u9- 9K9 (493) 

]^9'M9' = i(u, + K,j9j)Ml{uk + Kki9i) 
= ^-uM^u + ^{uiMlKki9i + Kij9jMluk) + ^Kij9jMlKki9i 
= ^uM^u + ^{uM^K9 + 9KM^u) + ^9kM^K9 
= ^uM^u + ^{uM'^K9 - uKM^9) + ^9KM^K9. 



^9'M^9' + 9'9 

= -uM^u + -u{M^K - M^K + 2)9 + 9K9 + -9kM^K9. (494) 



Here, we require M^K - M^K + 2 = 0. 
Using M = — M, we have 



which gives 



or 



2M''K=-2, 



K = -{M^)-'^, (495) 



^=(Mt)-i. (496) 



Here, we used K = —K, which is required for ii' as a new kernel. 
Finally, we have 

-9'M^9' + 9'9 = -uM^u + -9X9. 
2 2 2 

Therefore 

■ \ek0 



(497) 



= det(M)V2el^(M*)-^^, (498) 
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which is Eq.(226). Here, we used Eq.(179). 

Now, we continue to prove Eq.(228), where complex Grassmann variables are used. 
For complex variables. 



u = -^{ui + iu2), u'' = -^{ui - iu2), (499) 
V 2 v 2 



and 



ui = — ^(w + u^), U2 — —^{u — u^). (500) 
V 2 v2i 

The dual vector for ^(ui, U2) = {ui,U2\^) is given by 

*[wi,U2]= y"£'wiL'w2e"^"^e"^"^*^[ui,W2]. (501) 

Now, 

UiUi + U2U2 

= uu^ + u^u = u^u — u^u. (502) 



Thus, 



*[u,ut]= j DuDu^^u,u^)e^^''-'^^'' (503) 



which is Eq.(228). 

When *[u,ut] = e"^^" , ^^u,u^ = e"^^^" . 
Thus, we have 



which is Eq.(230). 



= Ce"^^", (504) 



C. Normal-Ordering of Hi 



A normal-ordered product of field operators is very convenient for calculation of the ground state energy (effective 
potential), because its expectation value vanishes identically Therefore before carrying out calculations on matrix 
elements, Eqs. (438)^^(442), it is necessary to normal order Hj. (j)^ can be rewritten using Eq.(428), 



where 



(27r)^/2 (2/(p))i/2 

ip + ^-+^p+, (505) 
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Using the above expressions, we readily obtain 

d^p 1 _ 1 
{2tt)D 2f{p) ~ 2 

(<Px - (fi) 





G3+1 

J 




- f 


{(Px - 








{4>x - 





ifi)^ : +6[ip-,ip+] : (<^x - 'P) ■■ +3[¥'", 



(507) 

(508) 

(509) 

(510) 
(511) 



Hi is now rewritten 



1 



+ 2 / loif^') 



+ (0x - 



An T 

D 



-m 



+ ^^(0c.-¥')^ + ^('/'x-^)^ 



(512) 



Substituting Eqs.(507)'~(511) into (512), we obtain 



Hr 



+^<^:(<^x-^)^:+^:(<^x-^)^:}, 



A 



(513) 



which is Eq.(446). 
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